THE GENERAL THEORY OF A CLASS OF LINEAR PARTIAL 
q-DIFFERENCE EQUATIONS* 


BY 


C. RAYMOND ADAMS 


INTRODUCTION 


During the past fifteen years a great deal of work has been done in 
building up a general theory of the linear difference equation in one in- 
dependent variable, or linear ordinary difference equation.t No study seems 
to have been made, however, of the linear difference equation in more than 
one independent variable, or linear partial difference equation. 

Essentially two distinct types of ordinary difference equations have been 
subjected to investigation: (i) that usually spoken of simply as the linear 
difference equation,t and (ii) that which has received the name of linear 
q-difference equation. These suggest three types of linear partial difference 
equations in two independent variables, which (if we confine our attention 
to homogeneous equations) may be written as follows: 


m n 


(4) (x,y) g(a+m—p,ytn—v) 


p=0 v=0 


m n 
(B) Any (x,y) 


v=0 


m 


n 
Any (2, y) y) 0, 


v=0 


in which gq and yr are constants, real or complex and different from zero; 
the a,,(#,y) are known analytic functions of x and y; and g(x,y) is the 

* Presented to the Society, December 28, 1923. 

7 For a comprehensive survey of what has been done, see Nérlund, Newere Untersuchungen 
iiber Differenzengleichungen, Encyklopadie der mathematischen Wissenschaften, vol. IIs, 
No. 6, Leipzig, Teubner, 1923, pp. 675-721. 

} Cf. Birkhoff, General theory of linear difference equations, these Transactions, 
vol. 12 (1911), pp. 243-284. 

§ Cf. Carmichael, The general theory of linear q-difference equations, American Journal 
of Mathematics, vol. 34 (1912), pp. 147-168. 
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function to be determined. To the equations (A), (B), and (C) we give 
respectively the names linear partial pure difference equation, linear partial 
q-difference equation, and linear partial difference equation of the intermediate 
type. If in (A), (B), and (C) mw be set equal to » and m equal to n, the 
following important classes of these partial difference equations are obtained: 


ay (x, 0, 


n 
> (x,y) x, y) = 0, 
v=0 


n 

av (x,y) y) = 
v=0 

It is the purpose of the present paper to develop a general theory of 
equations of the class (B’); equations of classes (A’) and (C’) will be 
treated later. For the most part the methods used are along the lines of 
those employed by Birkhoff* in his work on ordinary difference equations. 
In § 1, formal series solutions are found; § 2 is devoted to proving the 
convergence of these series. In § 3 are considered certain cases to which 
the methods of §$ 1, 2 do not apply. The object of § 4 is to indicate that 
when |qg| = |r| = 1, the equation in general has no analytic solution. The 
periodic functions which may be defined in terms of two solutions of an 
equation of this type are discussed in § 5. In §6 a theorem is proved 
which is in part an inverse of the results obtained in the earlier work. 

It should be observed that the chief results of the paper are immediately 
extensible to the case of an equation of the class (B’) in N independent 
variables; i. e., 


7 4 y 


v=0 


The author wishes to express his sincere thanks to Professor George 
D. Birkhoff, whose valuable suggestions have been of aid in the preparation 
of this paper. 

1. THE FORMAL SERIES SOLUTIONS 


The equation of the nth order of the class (B’) will be defined as 


n 
(1) > pr(x, y) g(qr-" x, = 0, 
vy=0 


* Loc. cit. 


n 
(4’) 
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in which g and y are real or complex constants not zero; the p,(2, y) 
are arbitrary known functions of the complex variables x and y and 
p(x, y) £0; and g(x,y) is the function to be determined. It is important 
to observe first of all that 

Every equation of the nth order of the class (B’) can be solved by deter- 
mining the solutions of an associated system of n equations of the first order 
of the class ( B’). 

In fact if we define 


g(x,y) = g(qi*a, ry) (i= 1,2,... 


the equation (1) is seen to be equivalent to the system 


ry) = 
(2) 
Pn-1 (az, y) 
n(2,y) go (x,y) 


_ Pala, y) 


gn ry) = po(x,y) 


Pr (x, y) 
— Yn . 
(x,y) 


We propose to show that an equation of the type (1) can in general be 
solved when the coefficients p, (x,y) are polynomials in x and y.* In that 
case the coefficients in the equivalent system (2), with which we shall deal, 
will be rational functions. It is unnecessary, however, for us to develop 
a theory for the solution of a system of equations with rational coefficients, for 

The solution of a system of n equations of the first order of the class (B’) 
with rational coefficients can in general be made to depend upon the solution 
of a second such system with polynomial coefficients. 

Proof. Let 


n 
gi(qz, ry) = y) gi(@,y) (i= 1,2,...,m), 
j= 


or in the matrix notation 
(3) G(qz,ry) = R(z,y) G(z,y), 


be the given system. With no loss of generality it may and will be 
assumed that each rational function rj(xz, y) is expressed as the quotient of 


* The problem is clearly the same as if the coefficients were taken to be rational functions. 


21* 


1), 
= 


286 Cc. R. ADAMS [July 
relatively prime polynomials. If d(z,y) is the least common denominator 


of the functions y¥), we have 


1 
R(z,y) = P(z,y), 


where P(x,y) is a matrix of polynomials. It will appear in the course 
of our later work (ef. §§ 2, 3) that, d(x, y) being a polynomial, solutions 
can in general be found for the equation 


g(qz,ry) = d(a,y)g(z,y). 


Let one such solution be denoted by ®(z,y) and set 


1 
G(z,y) O(z,y) G(x, y); 


then equation (3) becomes directly 
G(qz, ry) P(z,y) G(x,y). 


It is the chief purpose of the present paper to find analytic solutions 
of the system 


n 
(4) gi(qz,ry) = = g(x,y) 
J= 


where pi(x,y) is a polynomial 

py (2, y) = put pyr y + +++ + piitm zy”. 
In the matrix notation the system (4) is written 
(5) G(qz,ry) = P(z,y) G(z,y). 


Iteration of this equation suggests as symbolic solutions the two infinite 
products of matrices 


P(x, y) ry) r*y) 


Certain modifications of these will play an important part in our work. 
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Among all the polynomials py(z,y) there will be one (or more) whose 
degree in x is the greatest; let that degree be mw: similarly there will be 
one (or more) whose degree in y is the greatest; let that degree be »v. 
Then the following determinantal equations will be called the characteristic 
equations: 


for (0,0), | pij — Oy 0; 


for 00), | Pijpy — 94 


for (0, oo), | pijow — 


for (2, 0), | Pijpo — dy 


The roots of these equations will be designated respectively by 9;, 0, o/’, 
and o;” (4 = 1, 2,..., 
Our attention will at first be confined to the case in which the roots 


o;, 0, o;, and o;’’ are all different from zero.* We shall make the further 
restriction that 


(i+)), 


for any & and y which are positive integers or zero (but not both zero). 
Under these conditions and if we have |q|, |r|+1, and if the roots of 
each characteristic equation are distinct, it is found by direct substitution 
and reckoning of the coefficients that the system (4) is satisfied formally 
(a) by each of m sets of functions 


s1j (x, y) = [815 + 2+ sijor y + 


Snj (x,y) = + Snjor y + 


(j =1,2,...,m), 


*A method by which solutions can frequently (in fact “in general” for the case of 
a single equation) be obtained when this condition does not hold will be given in § 3. 


— 0: 
J, 
' 
oj gj 
. . . . . . . . 
(6) 
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where the g; are determined by the relations 


(b) by each of » sets of functions 


—t)/2 —Tt)/2 S1 ju $1j01 
(7) 


P—t)/2 Snj10 Snjol 
=1,2,... 
in which 
logg ’ 


and the @; are determined as 
_ logo. 


logq ’ 
(c) by each of nm sets of functions 


” 
” $1 j 


” 
(Gj = 1,2,...,%), 


where 
_ loge 
’ 


and (d) by each of m sets of functions 
si} (x,y) = whi | st + y+ 


(e—b Atd 8nj10 we 
Snj (x,y) q” xls + srjo y+: 


(j = 1,2,...,%), 


(July 
] 

qi =o or = = 
(8) 


1924] DIFFERENCE EQUATIONS 


in which 
yer 


log 
log q 


SJ 


Furthermore the sets of functions (6), and similarly (7), (8), and (9), are 
linearly independent. Throughout our work those determinations of log g 
and log r will be taken in which the coefficient of V —1 is positive or 
zero and less than 27.* 

We shall not assume, however, that the roots of each characteristic 
equation are distinct, but we shall make the hypothesis that the formal 
series solutions (6), (7), (8), and (9) do exist and that they are such that the 
following determinantal inequalities hold: 


= |ey| + 0, a’ = + 0, 
a’ = #0, a” = $0. 


The » sets of formal solutions (6) are then linearly independent and 
constitute a matrix S(z,y) which is a formal matrix solution of the 
equation (5). The same statement may be made for the solutions (7), (8), 
and (9); let the respective matrices be denoted by S’(z,y), S’’(z,y), 
and S’’’(a,y). It will be demonstrated (a) that if |qg| and |r| are both 
greater or both less than unity, the series in the matrix S(x,y) [S’ (zx, y)] 
converge uniformly and absolutely in the neighborhood of any place (2, y) 
for which x and y respectively lie inside [outside] certain associated circles 
72 Re] about the origins in their respective planes; and (b) that 
if one of the quantities |q| and |7| is greater and the other less than unity, 
the series in the matrix S’’(x,y) [S’’’(x,y)] converge uniformly and 
absolutely in the neighborhood of any place (x,y) for which x and y lie 
respectively inside [outside] and outside [inside] certain associated circles 
ft, | R,, rg] about the origins in their respective planes. 
The most general analytic matrix solution of (5) is given by 


H(z, y) G(x, y) A(z, 


*It is to be noted that the power of x occurring immediately before the bracketed series 
in each set of formal series solutions is not a uniquely determined factor; in place of x” 
in (6), for example, we may equally well use #5 or (ay)?s or, in fact, any function f(x,y) 
such that f(qx,ry) = cf (x,y), where c is a constant. 

7 We shall find it convenient to think of the complex variables x and y as represented 
by points in two distinct planes. In conformity with the usual convention a pair of values 
of w and y will be spoken of as the place (x,y). 
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where G(x,y) is any particular analytic matrix solution, and A(z, y) is 
any arbitrary matrix of analytic functions which possesses the multiplicative 
period system* (q, 7) and whose determinant | A(z, y)| is not identically zero, 


2. EXISTENCE THEOREMS 
Of the four cases |g| 2] 1, |r| ] 1 which are to be considered in this 


section, that in which both absolute values are greater than unity is wholly 
typical. We shall therefore confine ourselves in giving detailed discussion to 


The Case of \q\, |r| > 1. 


It has been observed that the matrix S(x,y) is a formal matrix solution 
of equation (5); that is, that 


S(qz,ry) = P(a,y) S(a,y). 


Now the determinant |S(z,y)| may be written 
+ dyoa + doy + 


The element in the ‘th row and jth column of the inverse matrix is the 
quotient of the cofactor of the element in the jth row and 7th column of 
this determinant by the determinant itself, and is therefore given by a series 


(11) y) == Pi siz + Sino + Sijor y+ 


Let T(x,y) denote the matrix obtained by breaking off the elements 
of S(x,y) so as to retain in the series only the terms of degree less than 
or equal to k—1, or more generally by replacing sj(z,y) by ty(z, y), 
where the series in the latter is convergent for (2, y) in the vicinity of 
the place (0,0) and has the same terms as the series in syj(x,y) up to and 
including those of degree k—1. We may then define a matrix Q(z, y) 
by the relation 


T(qx,ry) = Q(a,y) T(2,y). 


* The matrix A(z,y) shall be said to possess the multiplicative period system (q,r) 
if A(qz,ry) = A(a,y). 
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The matrix Q(x, y) is a matrix of functions gjj(2, y) each of which has an 
expansion in powers of x and y that agrees precisely with p(x, y) up to 
and including terms of degree k—1. This follows upon comparing 


= T(qz,ry) T (a, y) and P(x,y) = S(qz,ry) S'(2, y); 


for T(qx,ry) is the same as S(qz,ry) up to and including terms of 
degree k—1, and similarly for 7’—'(x, y) and because ty (a, y) 
of T(x, y) is given by (11) up to and including terms of degree k—1. 
We therefore have 


P(x,y) = Q(2,y)+ M(z,y), 


where M(x, vy) is a matrix of power series in x and y whose lowest degree 
terms are of degree k or higher and which converge in the vicinity of 
the place (0,0). Then N(2z,y), defined by 


M(x2,y) = Q(2,y) N(a,y), 
is a matrix of power series of the same type as those in M/(z,y). Hence 
P(x,y) = Q(2,y) N(2,y)], 


in which J is the unit matrix and N(x, y) is a matrix of power series each 
of which contains no term of degree less than k and which converges in 
the vicinity of (0, 0). 

THEOREM I, A. Form the product of matrices 


> 


Each element of Pm(x,y) converges, for k sufficiently large, to a definite 
limit function uj (x,y), independent of k, as m becomes infinite. This function 
is analytic throughout the entire finite x- and y-planes except perhaps at 
«= 0, and in this region is identical with the corresponding element sy(x, y) 
of S(x,y). 

Proof. We may write 


Pm (x, y) T(x,y) Pn(x,y), 
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where 


|r (=, 
q 


The elements of 7'(x, y) are polynomials or series convergent in the vicinity 
of (0,0) (multiplied by a power of x). In order to show that the elements 
of Pm(x,y) converge, it will be sufficient to prove that those of Pm(a,y) do. 
We proceed to a proof of this fact. 

The matrix P,,(x,y) is the product of matrices of the type 


P(x,y) = I14+T (a2, y) N(x, y) T(z, y) 
The second term is a matrix 
(a “da (a, y)) = (04 (x, y)) = O(a, y) 
where the 4; (x,y) stand for power series in x and y whose lowest degree 
terms are of degree k or higher and which converge in the neighborhood 


of (0,0). 
The matrix Pm(x2,y) may now be written 


The 7th element in the jth column (7,7 = 1,2,...,m) is 


m 

x 

Pm:ij (2, y) — + (= 


+E 


+1 q r 


6;; denoting the element in the th row and jth column of the unit matrix. 


q”" ym 
~ L y | y 
q a ? “a 
y 
yh 
(12) 
m—t 
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Let it next be noted that the product of two (and hence of N) matrices 
of the type of O(2,y) which differ only in the series 4 is another matrix 
of the same type. For, if the two factor matrices are O(z, y) and O'(z, y) 
and their product is 9” (a, y), we have (denoting by [4] a series of the type 2) 


+ Biz 05; + + Onj 


p 


™ 


This remark makes it clear that we have 


m 


(13) 
m—1 


m r 


If p is equal to the smaller of the two quantities |q| and |r| and if 2 and y 
are inside certain associated circles r,, 72 about the origins in their respective 
planes, then the typical element in the /th term of the right-hand member 
of (13) is less in absolute value than 


M! 


because when (x,y) is in this region all the functions 4(x,y) are less 
than some constant M. Furthermore if s is a positive integer assigned 
arbitrarily, we can take & so large that the inequality 


> pe 
will hold for «, 8 = 1,2,...,. Let any positive integer s be assigned 


and a suitable k be found; it follows that the series (13) will then be less 
term by term in absolute value than the series 
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As m is allowed to increase without limit this series approaches the value 


(1 


The ratio 


since p is >1 and sis a positive integer, is less than unity. The element in 
the 7th row and jth column of Pm (x,y), after being multiplied by afi? 4, 
therefore converges uniformly and absolutely to a function analytic in the 
neighborhood of places (x,y) such that |z| is <r, and |y|<~rz. Hence 
each element of Pn(x,y) converges in a similar manner in the neighbor- 
hood of places (x, y) such that 0+ |y| <7. 

We may write 


P(x, y) P| 


and apply the argument as above to Pm—«(x/q', y/r*). It follows, since |q| 
and |r| are >1, that if any place (x,y) be assigned where x and y are 
in the finite parts of their respective planes and x +0, ¢ can be chosen 
so large that (a/gt) (+0) and y/r* will lie within the respective circles 7, 
and 72. We conclude that the elements of P»(2,y) converge to functions 
uij(x,y) analytic over the entire finite parts of the z- and y-planes except 
at z= 0. 

The functions w(x, y) are independent of k (so long as k is sufficiently 
large), because if k’ is any second value for k and 7"(xz,y) the matrix 
corresponding to 7'(x,y), we may write 


r 


Pras y) = 4) 


ps pons p’ 
4) Pat (5,4 
Y\ YI... 
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Then if the elements be expanded in a sum like (12) and m be allowed 
to become infinite, the resulting multiple series for uy(a2,y) will be term 
by icrm identical with that obtained above. 

It remains for us to show that if x and y are anywhere in the finite 
parts of their respective planes and x +0, then u(x, y) = sy(x,y). We 
shall first show this to be the case if x(+0) and y are respectively within 
the associated circles 7,, v2. Let (x,y) be any place in this region. For 
this pair of values the limit function %j(a,y) of the 7th element in the jth 
column of Pm(a,y) differs from dj by a quantity less in absolute value than 


n p 


which is itself less than 


nM 


since 1/p* is <1. Hence we have 
| y) — < 


If s (and with it k) be allowed to increase without limit, it is clear that 
y) approaches Therefore is 


uy (x,y) (a, y) (i,7 = 1, n), 


for we have shown wj(x,y) to be independent of k, and as k becomes 
infinite ty(x, y) becomes syj(xz, y). Thus the function w;(x, y) is 
represented by sj (a, y) for + |a|<r, and |y|<r,. But the function 
uj(x,y) is analytic without singularities of any kind over the entire 
finite a- and y-planes except at « = 0. Hence the power series in s; (x, y) 
converges for all finite values of x and y and represents x "i wy (a, y) 
throughout this region; it follows that wy(z, y) = sy(x, y) in the 
region specified. 
The proof of the theorem as stated is thus complete. 


But 

= ? 

| | 
n 
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We may consider in a similar manner the matrix S’ (x,y). It hag 


already been noted that this matrix is a second formal solution of the 
matrix equation (5); that is, that 


S'(qz, ry) = P(x,y) S(2,y), 
which, for our immediate purposes, may be written more conveniently as 
S(x,y) = P'(a2,y) S ry). 
The determinant | S’(z, y)| is of the form 


, 


nu(t?—t/2 nv(r?—r)/2 p{+p$+---+! dio do. 

y 


The element in the 7th row and jth column of the inverse matrix 


consequently seen to be 


Y 


If 7" (x, y) denotes the matrix obtained by replacing sj; (x, y) in S’ (x, y) 
by tj; (x, y) where the series in the latter is known to be convergent in 
the neighborhood of (oo, co) and has the same terms as the series in 
sij (a, y) up to and including those of degree k—1 in 1/x and 1/y, then 
we may define a matrix Q’—! (xz, y) by the relation 


T' (x,y) = y) T' (qz, ry). 
By comparison of 


Q' (x,y) = T(z, y) T’ (qa, ry) and P(x, y) = (x2, y) (qz, ry), 


we see that Q’~!(z, y) is a matrix of functions gj; (x, y) of the form 
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in which the series has exactly the same terms of degree less than / in 
and 1/y as has y” pY¥ (x,y), (a, y) standing for the element in 
the ‘th row and jth column of the matrix P-'(z,y). It follows that 


P (2, y) = (2, (2, y)], 


where the elements of the matrix N’ (x, y) are power series in 1/x and 1/y 
which contain no terms of degree lower than k in these variables and which 
are convergent in the vicinity of (oo, o). 

We may now state 

THEOREM I, B. Form the product of matrices 


I, (x, y) = P— (x,y) ry)--- 2, y) r™y). 


Each element of Im (x, y) converges, for k sufficiently large, to a definite 
limit function v(x, y), independent of k, as m becomes infinite. This 
function is analytic throughout the finite x-and y-planes with the following 
exceptions: 2 = 0, y = 0, and places which are poles of an element of 
one of the matrices y), P(qaz,ry),---. For x and y outside 
certain associated circles about the origins in their respective planes and 
x,y + ©, vy(x, y) is identical with the corresponding element sijj(x, y) 
of S' (a, y). 

The proof of this theorem runs entirely parallel to that of Theorem I, A; 
we shall omit the details and pass on to the statement of the corresponding 
pair of theorems which hold in 

The Case of |\r|<1. 


Here the two symbolic solutions of the equation (5), P(z/q, y/r) 
‘P(z/q?, y/r*)--- and P— (a, y) P— (qx, ry)---, interchange the réles which 
they played in the discussion of the preceding case. The theorems follow. 

THEOREM I, C. Form the product of matrices 


In(x,y) = P(x, y) (qa, ry)--- P— (g™— 2, r™—1y) T(q™z, rmy), 


Each element of Hm(x, y) converges, for k sufficiently large, to a definite 
limit function uj(x, y), independent of k, as m becomes infinite. This 
function is analytic throughout the finite x- and y-planes except perhaps 
at x = 0, and except at places which are poles of an element of one of 
the matrices P—(x, y), P—(qx, ry),--+. For x and y inside certain 
associated circles 11, 12 about the origins in their respective planes and 
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x +0, w(x, y) a identical with the corresponding element sy(x, y) of 
S(x,y). 
THEOREM I, D. Form the product of matrices 


Each element of Pm(ax, y) converges, for k sufficiently large, to a definite 
limit function v(x, y), independent of k, as m becomes infinite. This 
Junction is analytic throughout the finite x- and y-planes except at x = 0 
and y = 0, and in this region is identical with the corresponding element 
s(x, y) of (ax, y). 

Of course the wy(a2, y) and vy(, y) of this case are entirely distinct 
from the solutions of the preceding case for which the same notation 
was used. 


The Cases of |\q|>1, |r| <1 and |q| <1, |r| >1. 

It is in the discussion of these cases that the matrices S”(a, y) and 
S’’ (a, y) are of importance. We define a matrix (x, y) [T’”’ (2, y)] 
by replacing the series there by power series in x and 1/y [1/a and y] 
which are known to converge in the vicinity of the place (0, 0) [(o, 0)] 
and which agree precisely with the series in S” (zx, y) [S’’ (x, y)] up to 


and including terms of degree k—1 in x and 1/y|[1/z andy]. Proceeding in 
entirely the same way as before we can prove, in each of the two cases, 
two theorems analogous to those above. We shall not pause to state the 
theorems here; it will suffice to remark that in both cases the series in 
the matrix 8’ (x, y) [S’’ (a, y)] converge for x and y respectively inside 
[outside] and outside [inside] certain associated circles 7r,,R, [R,, rz] about 
the origins in their respective planes, and that by means of the equation (5) 
itself, analytic or meromorphic solutions are defined throughout a more 
extended region. 

In the existence theorems for each of the four cases, use is made of 
but two of the four formal matrix solutions; for example, if |q| and |r 
are both greater than unity, we have used only S(z, y) and S’(z, y). 
It would be of interest to know just what is the significance of the two 
remaining formal solutions, 


3. THE LINEAR EQUATION OF THE FIRST ORDER WITH POLYNOMIAL COEFFICIENT 
FOR THE CASE IN WHICH A CHARACTERISTIC EQUATION HAS THE ROOT ZERO 
It will be recalled that in § 1 one of the restrictions which it was 

found necessary to make in order to obtain the formal solutions (6), (7), 


Pm (a ’ y) I ’ Fa I g’ 1 ~ 
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(s), and (9) is that the roots of each characteristic equation be all different 
from zero. An important particular case is that in which the system (4) 
consists of a single equation (i.e., when » = 1): 


g(qz,ry) = plz, y) g(a, y), 


p(x, y) = pt pw s+ ety”. 


In this case if p = 0, the one root of the characteristic equation for 
(0,0) vanishes and no formal series solutions of the type (6) can be found. 
Then if |g| and |7| are both > or <1, one of the two solutions that are 
usually obtainable is lacking. We desire to show, however, that a solution 
ean in general be found to take its place. We shall restrict our discusssion 
entirely to this case, although the method introduced is of rather general 
application. It can be used directly if |g| and |v! are both > or <1 
and a solution is missing because p,, = 0, and with slight modification 
if q| is >1 and |r| is <1 or |q| is <1 and |r| is >1 and a solution 
is missing because py, = 0 Or pyo = O. It may also be possible at times 
to use the method effectively when n is >1 and one (or more) of the 
roots of a characteristic equation vanishes. 


We shall accordingly assume |q!, |r|>1 or |q|, |r|<1 and proceed 
to obtain a solution for the equation (14) when p = 0. If p(x, y) can 
be factored into aX yf (where x and 4 are positive integers or zero, but 
not both zero) multiplied by a polynomial p’(«, y) whose constant term 
is not zero, then a solution of (14) can be obtained by solving the two 


equations 


= (ay), = p'(2,y) g(x,y), 


and multiplying the solutions together. The first of these equations is 
satisfied by 


logr _ logy 
logqg’ logr’* 


y) = where t = 


As for the second equation, that is solvable at once by the methods of 

A geometrical interpretation of what it means for p(z,y) to be factor- 
able into «*y* p'(x, y) is worthy of notice. The terms of the polynomial 
p(x, y) are of the form «”y”; let the point corresponding to each term 
be plotted in the m, m plane. Clearly each such point will lie in the first 
quadrant or upon the positive m- or n-axis, and a necessary and sufficient 


9 


(14) 
| 
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condition that p(az, y) be so factorable is that there be one point of this 
set in the m,n plane such that if there be drawn through it parallels to 
the axes, every other point of the set will lie above or upon the paralle| 
to the m-axis and to the right or upon the parallel to the n-axis. 

When p(az, y) cannot thus be factored, our method of solving equation (14) 
will be to effect upon the variables 2 and y a transformation which wil] 
leave the given equation invariant in form but will change p(z, y) intoa 
polynomial p(x, y) which can be factored in the desired manner. We 
propose to show that 


(15) 


where each of the exponents is a suitably determined positive integer or 
zero, is such a transformation. It is readily seen that under this trans- 
formation the equation (14) becomes 


(16) 


in which 


Hence it is our purpose to prove the following 

THEOREM. The quantities a, B, y, 6 in (15), each of which is a positive 
integer or zero, can in general be so chosen that 

(a) p(x, will be factorable into x* y* p' (z, y), where (x, y) has 
a constant term different from zero, and 

(b) |q| and |r| will be both greater or both less than unity. 

Proof. The effect of the transformation (15) on the polynomial p (z, y) 
may best be observed in the m,n plane. By (15) the term 2” y” goes 
into the term 


zm y” — Zamtyn | 


Thus the transformation of the m,n plane which corresponds to (15) is 


m= am+yn, 


Bm + én, 


J = ad—By $0. 


2 = 
| ad—By + 0, 
y= ay’, 
| ry) = P(z, y) g(@, Y), 
a 
q’ 
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We shall think of this as a transformation of the axes, leaving the points 
of the plane fixed. Then the equations, referred to the m, n axes, of 
the im, axes are respectively 


These lines lie in the second and fourth quadrants, or at most, one of them 
may coincide with the m- or the n-axis. Furthermore, since for m and n 
positive m and % are also positive, it follows that the angle @, less than 7, 
between the positive m- and y-axes must include as a part of itself the 
angle between the positive m- and n-axes. By a proper choice of the 
exponents in (15) the m- and 7i-axes can be put into any desired position 
of this type so long as the slopes of the m, m axes are rational (or in 
particular, infinite). 

Let us now consider again our set of points in the m, m plane corres- 
ponding to terms of p(x, y). A convex polygon can be drawn such that 
each vertex is a point of the set and every point of the set lies within 
or on the boundary of the polygon; : 
let this polygon be constructed. ‘ 
Clearly it will have a vertex whose 
m-coordinate is less than that of any 
other vertex, or of any other but one. 
Call this point A; in case there are 
two such points, call the one with 
the larger n-coérdinate A. Starting 
from this point proceed around 
the polygon in a counter-clockwise 
direction, lettering the successive 
vertices B, C,... until a vertex is 
reached whose n-coérdinate is less 
than that of any other, or of any 
other except one. Let this last vertex be called N; or, if there are two such 
points, call the second one met (the one with the larger m-coérdinate) N. 

Now the slopes of all the lines AB, BC, ... are rational and negative, 
or in particular infinite or zero. If the slope of AB is not infinite, draw 
through A a parallel AR to the n-axis; if the slope of MN is not zero, 
draw through N a parallel NS to the m-axis. Then it is possible — 
given any vertex of the set A, B, ..., N (exclusive of A if AB has an 
infinite slope and exclusive of N if MN has a zero slope), say B — to 
choose @, 8, y, 6 so that the 7- and m-axes [the m- and 7-axes] will 


| B a 
n j m, = m. 
| 
= 
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be parallel respectively to AR and AB, or to BA and BC, or..., or 
to NM and NS. Any of these choices will yield a transformation (15) 
which will take p(z, y) into a polynomial of the desired factorable type, 

Secondly consider the condition (b). By hypothesis we have | q|,|7| > 1 
or |q|,|r|<1, and it is necessary that J be positive or negative. Let 
us assume |q|,|7| > 1 and try to determine a, 8, y, 0 so that J will be 
> 0 and condition (>), as well as condition (a), will be satisfied. It is 
readily seen that |q/| and |r| will both be greater [less] than unity if 
and only if 


a _ log\g log|q| 
log |r| 


log|r| 


The second of these possibilities, however, necessitates J being < 0. Hence 
if |g| and |r| are >1 and we are to use a transformation (15) in which 
J is >0O, a necessary and sufficient condition that the transformation 
satisfy the requirement (b) is that @, 8, y, 0 be determined to satisfy 
the inequality 


_ log\q| _ B 
y~ logir|~ 6 


But log|q|-+log|r| is a positive constant. Its negative therefore lies, in 
general, between the values of the slopes of two consecutive lines of the 
set RA, AB, BC,..., MN, NS; in particular, its negative may be equal 
to the slope of one of the lines AB, BC,..., MN (AB excluded if its 
slope is infinite and MN if its slope is zero). In the general case, if 
a, 8, vy, 8 be so chosen as to make the 7- and m-axes parallel respectively 
to the two lines of the set RA,..'., NS between the values of whose 
slopes — log|q| + log|7| lies, then the condition (b) as well as the condition (@) 
will be satisfied. In the particular case a solution cannot be obtained by 
this method. 

Several observations should be made at this point. First, by determining 
a, 8, y, 0 so that the réles of the m -and 7-axes in the preceding paragraph 
are interchanged, we obtain a transformation (15) satisfying (a) and ()) 
in which J is negative and | q| and |7| are both less than unity. Secondly, 
in the case of |q|, |x|<1 a wholly similar discussion may be given. 
Thirdly, it is clear that in either case we can make | q|,|7| > 1, so that 
the words “or both less” may be stricken from part (b) in the statement 
of the theorem. 
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A solution for the equation (16) can be found by solving the two equations 
TY) = gar, ry) = y) y) 


and multiplying the solutions together; such a solution will be analytic 
throughout the finite z- and y-planes except at x = Oand y = 0. Trans- 
forming this solution into terms of 2 and y by the inverse of trans- 
formation (15) we obtain a solution of (14) which is analytic throughout 
the finite z- and y-planes away from x = 0 and y = 0. 


4, THE CASE oF |q| = |r| = 1 

THEOREM. When |q| = |r| = 1 the equation of the class ( B’) in general 
admits no analytic solution that is not identically zero. 

We distinguish two sub-cases: (i) when at least one of the constants q, r 
is not an mth root of unity for any positive integer n; (ii) when both q 
and » are roots of unity. The proofs will not be given here, for in each 
case a demonstration can be effected by the methods used by Carmichael* 
in proving the corresponding theorem for the linear ordinary q-difference 
equation; no very essential modification of Carmichael’s work is needed 
beyond the introduction of a second independent variable.+ 


* Loc. cit., pp. 165-168. 

+ One point alone deserves special mention. In constructing a proof in subcase (i) it 
becomes necessary to show that, any positive integer M being given and any 2 and y 
being assigned, an integer » can be found which is > M and such that the place 
(q*z,r"y) is as near as you please to the place (x,y). This problem is equivalent to 
that of finding an » > M such that we shall have 


and |r*>—1|<e, 


¢(>0) being preassigned and arbitrarily small. If we set « = argg and 8 = argr, 
this reduces to determining integers n (> M),m, and p to satisfy at the same time the 
two inequalities 
|na—m-22|<e' and 

That such a determination is possible follows immediately from a generalized form of 
Theorem III, pp. 7-8, of Minkowski’s Diophantische Approximationen, Leipzig, Teubner, 1907. 
If in Minkowski's proof (p. 7) we allow his variable z to take on the values0, M,2M,...,@M, 
where M is an arbitrarily assigned positive integer, the rest of the proof goes through 
without change, and we obtain a theorem which may be stated as follows (for M = 1 
it is Minkowski’s): Given two arbitrary real numbers a and b and any two positive 
integers M and t; there can always be found at least one integer n, M<n=<#M, and 
two other integers m and p such that the following inequalities hold simultaneously: 


ina—m| <4 
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5. THE PERIODIC FUNCTIONS 

The discussion of this section will be confined to the case of |q|,|r| > 1: 
each of the other cases covered by the existence theorems of § 2 may be 
treated in the same manner. 

It has been shown in § 2 that, under suitable restrictions, the matrix 
equation (5) admits of two matrix solutions, which may here be denoted 
for convenience by Soo(x, y) and S,,..(%, y) respectively. These solutions 
have been found to possess the following properties. Soo(z, y) is analytic 
when x (+ 0) is inside an arbitrarily large circle R, about the origin in 
the z-plane and y is inside an arbitrarily large circle R, about the origin 
in the y-plane, and is identically equal to S(x,y). In a similar manner 
S.co(x,y) is analytic save for poles when x (40) and y(+0) are 
respectively outside arbitrarily small circles r,, rz, about the origins in their 
respective planes. If (x, y) lies within the region of convergence of the 
formal series solutions (7), we have S,,..(x, y) = (a, y); if (a, y): lies 
outside that region, S.,..(2, y) is given by 


y) P(qa, ry) (qa, (qx, ry), 


where s is a suitable positive integer. 
We now define a matrix A(xz,y) by the equation 


(17) Soo(x,y) = A(z, y), 


from which it is clear that A(z, y) is a matrix of multiplicatively periodic 
functions such that 


A(qu,ry) = A(az,y). 


The elements of A(x,y) are thus analytic, except perhaps for poles, in 
the region consisting of the x-plane between the circles 7, and AR, and 
of the y-plane between the circles 7, and A; it will first be proved that 
no such poles can occur. 

We have A(z, y) = SZ}, (x, y) Soo(a, y). The elements of S51, (2, y) 
may have singularities at the poles of elements of S,,..(x, y) and at the 
zeros of the determinant | S,,..(z, y)|. But outside of certain definite as- 
sociated circles R,, Rz the elements of S,, (x, y), so long as neither z 
nor y is actually infinite, have no poles, and outside certain other definite 
associated circles R,, A, the determinant |S,,,.(z, y)|, so long as neither 
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, nor y is actually infinite, does not vanish.* Now let any place (2, y ) 
be chosen such that z, lies between 7, and R, and y, lies between r, 
and Rs. We have 


= Alan, ry) = = Alin, ry), 


and s can be taken so large that g’z, will lie outside the larger of the 
circles R,, A, and yet inside R,, and r*y, will lie outside the larger of 
the circles R,, A, but inside R,. Hence A(a, y) must be analytic at 
(q’a,, y1) and therefore at (2,,¥,). From this we conclude 

THEOREM II, AB. The elements of the matrix A(x, y) defined by (17) are 
functions possessing the multiplicative period system (q, 1); these functions 
are analytic in the region consisting of the x-plane between circles r, and 
R, about the origin and of the y-plane between circles rz and Ry about the 
origin, where 1, and rg are arbitrarily small and R, and Rz are arbitrarily 
large. 

We shall find it worth while to examine these periodic functions a little 
more in detail. In order that they may be single-valued functions of position 
it will be convenient to think of each of the variables x and y as free to 
range over an infinitely many-sheeted Riemann surface of the type used 
for the logarithm, with branch points of infinite order at the origin and 
at infinity. Let us then make the transformation 


which maps the Riemann surfaces for x and y respectively in a one-to-one 
and conformal manner on the ¢-plane and t-plane, and define 


A(x,y) = B(t,r). 


B(t,c) is thus a matrix of functions single-valued and analytic over the 
region made up of the portions of the ¢-plane and t-plane within arbitrarily 
large circles about the respective origins. Moreover the elements of B(t, r) 
have additive, rather than multiplicative, period systems. 

Let x now be made to traverse a positive circuit about the origin in 
its own plane (y remaining fixed). Syo(x, y) becomes Soo(x, y)K, where 


K stands for the matrix (2”— dy), 64 being as before the Kronecker 0. 


* Cf. inequalities (10). 


t 
x= y=r, 
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At the same time S,,,,(7,y) changes to 


where L is the matrix (ei 6,). Thus when x has made a positive 
circuit of the origin, A(z, y) is replaced by 


on /— og2 1 


If, then, aj(a,y) is the element in the ‘th row and jth column of A (2, y), 
the passage of x around such a circuit changes a(x, ¥) to 


27 1 yt 27% 1(- 


Similarly if y be caused to traverse a circuit about the origin in the 
positive sense (a remaining fixed), ajj(a,y) is replaced by 


y —onV—lvr on? 


It follows that if b(¢,7) denotes the element in the zth row and jth column 
of B(t,r), then 


bj(t+1,¢+1) = by(t,s), 


(18) by (t+ — 


— (— —1( logq bi ( ) ; 


— (—1 e b(t, t). 


Thus the elements of the matrix B(t,7z) are seen to be triply periodic, 
or perhaps one might better say “triply quasi-periodic.” As such these 
functions are a special case of the functions treated by Cousin* in a memoir 
in which he discusses, with only slight restrictions, the properties of the 
general meromorphic function of two complex variables whose zeros admit 


* Sur les fonctions triplement périodiques de deux variables, Acta Mathematica, 
vol. 33 (1910), pp. 105-232. 


[July 
by (t, 2nV—1 
log r 
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three systems of periods.* These restrictions are in general satisfied by 
the functions bj(t,c). The functions bj(t,7) are a more restricted type 
than those dealt with by Cousin in that they are entire, rather than 
meromorphic functions and in that the exponents of e in the equations (18) 
are linear functions of ¢ and t rather than the wholly general entire functions 
which appear in the corresponding equations of Cousin. The Cousin paper 
also shows how a triply periodic function can be expressed in explicit form 
in terms of series of the type 


+00 


in which ®,(¢) is essentially an exponential function of ¢ multiplied by 
a product of 6-functions of a (a” being a constant (7 = 1,2,..., 
Hence the functions );;(¢,7) are expressible in such a manner as this, which 
is remotely analogous to the explicit form found by Birkhoffv for the periodic 
functions that arise in the theory of the linear ordinary q-difference equation. 


6. A MORE GENERAL PROBLEM; AN INVERSE THEOREM 


In the four matrices of formal series solutions as they were found in § 1, 
there seems to be a lack of symmetry. This lack, however, is only apparent, 
for in defining our characteristic equations we might equally well have 
proceeded as follows. Among all the polynomials pi(z,v) there will be 
one (or more) which contains x[y] to a maximum power; let that power 
be «[v]; similarly there will be one (or more) which contains x [y] to 
a minimum power; let that power be *[4]. Then the characteristic 
equations are 


for (0,0), Pijxa — 94 


for (2,00), | — dj | 


for (0, oo), | o” 


yer 
| 


for (2, 0), | 


*The zeros of f(t,t) are said to admit the system of periods (a,b) if f(t+a,r+ 5) 

F(t, 2), where g(t,t) is an entire function. 

7 The generalized Riemann problem for ..., Proceedings of the American 
Academy of Arts and Sciences, vol. 49 (1913), pp. 561-564. 


| — 
| = 
— 
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and the formal matrix solutions are found to be 


y) 2 yA? tT)/2 [ sij Sij10 r + Sijor Y + 1), 


” 
y 


S(2,y) (a ? E + Sijor 


In this form there is evidently complete symmetry, and we note that 
when x = 4 = 0 the matrices reduce to those of § 1. 

It immediately becomes clear, after the preceding remark, that the 
polynomials py(a,y) may just as well be taken as polynomials in 2, 1/z, 
y, and 1/y. In that event the formal work of the foregoing paragraph is 
unchanged, but « and 2 are negative (i. e., if terms in 1/x and 1/y actually 
occur). Furthermore it should be observed that in the cases of |q|, |7|/>1 
and of |q|, |r| <1, existence theorems can be proved almost entirely as 
above if the elements of the matrix are any functions whatever which are 
expressible in each of the following forms, the first valid in the vicinity 
of (0,0) and the second valid in the vicinity of (00,0): 


[convergent power series in x and y], 


1 
a” y’ |convergent power series in — and 
x y 


(provided, of course, the restrictive conditions imposed in § 1 are satisfied); 
the only difference in the results is that the solutions obtained may be valid 
only in definitely limited regions about (0,0) and (a, o) respectively. 
It is clear that a similar remark applies in the two remaining cases, 
q| > 1, |r| <1 and |q| <1, |r| >1. 

We are now in a position to prove a theorem which is in part an inverse 
of some of the results that have been found in this and preceding 
sections. 


uly 
)) 
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III, ABCD. Let G(x,y) and H(ax,y) be two matrices whose 
cements, after multiplication by o-Ps in the case of 
Gir.y) and by q 22 in the case of H(x,y), are single- 
ralued functions analytic save for poles over the extended x- and y-planes 
with the exception of the places (0,0) and (a#,«@) and such that 


lim g(x,y) —A = sy, 
xr=0 

y=0 

lim q 


2 


? 


y=o 


where the s’s and @’s are constants subject to the conditions 


logr log r 


for any Band y that are positive or negative integers or zero (i+ 7); and 


where z, 2, and v are integers (or zero) and 


log y 
log 


Also let the two matrices be connected by the relation 


G (x, y) y) A(z, 
in which A(x, y) is a matrix of functions possessing the multiplicative period 
system (q,r). Then G(x,y) and H(x,y) are solutions of a system (5) in 
which the elements of P(x,y) are rational functions. 
Proof. By hypothesis we have 
H-(qx,ry) G(qz,ry) = H(a,y) G(2,y). 


Let P(x, y) denote the product of matrices 


G(qz,ry) (x,y) = H(qz,ry) H(2,y); 


qi, |rj>1 

log q 
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its elements are single-valued and analytic except for poles over the 
extended z- and y-planes apart from the places (9,0) and (@,o), 

Upon calculating the elements gij(x,y) [hy(x, y)] of @ (a, y) 
it becomes clear, by virtue of our hypotheses, that 


lim gy(x,y) oP = sij 
y=0 


in hi (x, y) q 7 


1 


where sj; [sj] is the element in the 7th row and jth column of the inverse 
matrix of (sj) [(sij)]. Hence, if (7, y) is in the neighborhood of (0, 0), 
we have 


y) = ry) y) = y/ (sic + Nije), 
« « 


where jije(i,j,@ = 1, 2,...,”) approaches zero as (xz, y) approaches 
(0,0). Similarly, if (2, y) is in the neighborhood of (oc, 0), we have 


n 
py 2, y) “ey jp q * ( Sic Saj + Nija) 
a=1 


in which yije(7,j,¢ = 1,3 n) approaches zero as (x, y) approaches 
(co, 0). It follows that the elements of P(x, y) are analytic or have 
poles at (0,0) and at (co, co). Consequently they are analytic except 
for poles over the extended z- and y-planes without exception, and are 
rational functions of x and y. 

We have, then, the two equations 


(20) G(qz,ry) = Gla, y), 
(21) H(qz,ry) = P(a«,y) H(2,y). 


From the preceding paragraph it follows that the elements of P(z, y) 
have at (0,0) and at (oo, co) the requisite forms (19). It remains to 
show that the roots of the characteristic equations for (0, 0) and for (oo, ©) 
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satisfy the conditions imposed in § 1. To do this let us rewrite (20) in 
the form 


n 


For any particular 7 this system provides n linear non-homogeneous equations 
for the determination of the functions y) = 
If they be solved for piu(x, y) and (x, y) be allowed to approach (0, 0), 
the coefficient in the leading term of the expansion for pi,(a, y) at (0, 0) 
being denoted by piuxza, we have 


= 


where |sjm?| stands for the determinant obtained from | sym| by replacing 
the elements of the uth row by q“sy(j = 1,2,...,m). (It is clearly 
possible thus to solve, because by hypothesis | sz, || 0.) The characteristic 
equation for (0, 0) is 


Piuzd— Sin | = 0; 


upon examination its roots are seen to be qg”, qg”2,..., q°». In the same 
manner we may show from (21) that, if pix... is the coefficient in the leading 
term of the expansion for pju(a2, y) at (oc, oo), then we have 


iupey 


and that the roots of the characteristic equation for (co, «), 


| ae 
| Piupy | 


are q?t, gft,..., g’». It then follows from our hypotheses that the matrix 
satisfies all the conditions that it has been found necessary, in our earlier 
work, to impose upon the matrix of coefficients in the system (5). This 
completes the proof of the theorem. 


a 
(iu) 
— 
Slm 
_ Slm 
| Slm | 
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If one of the two numbers |q| and |r| is greater and the other less 
than unity, an entirely similar theorem can be stated and proved. 

In conclusion, we may call attention to the fact that a discussion of 
the equation (2’) such as has been given covers also the cases of equations 
(A’) and (C”) for certain types of coefficient functions. In the first place, 
subjecting equation (5) to the transformation x qe ¥ r*, we see at 
once that our discussion disposes of the case of a system of equations of 
the first order of the class (A’) when the coefficients are polynomials (or 
such other functions as have been noted above) in q” and 7%. Secondly, 
if in (5) we seta = t, y r*|, we find that our work 
has covered the case of a system of equations of the first order of the 
class (C") with coefficients which are polynomials (or such other functions 
as we have noted above) in g® and y [a and r¥]. 

Brown UNIVERsITY, 
ProvipEncer, R. I. 


THE SUMMABILITY OF THE TRIPLE FOURIER SERIES AT 
POINTS OF DISCONTINUITY OF THE FUNCTION DEVELOPED* 


BY 


BESS M. EVERSULL 


In a previous paper? the author has studied the summability of the 
development of a function of three variables in a triple Fourier series at 
points of continuity of the function developed, using the method originated 
by Fejér and applied by him to problems involving simple series. For the 
latter series the proof of the summability at points of discontinuity of the 
first kind (finite jumps) is quite similar to that for points of continuity, 
and the two eases can be treated by means of a single discussion. In 
passing to the case of the double Fourier series, it is found that the 
study of the behavior of the series at points of discontinuity of an analogous 
type presents difficulties and complications that do not arise in connection 
with points of continuity.{ When we go on to the case of triple series, 
we find that the difficulties and complications of the corresponding problem 
are still further increased. 

In the present paper a study is made of the summability of the triple 
Fourier series at points of discontinuity of the type that would be apt to 
arise in physical applications. This includes discontinuities lying on plane 
or curved surfaces, and such that the function to be developed approaches 
the same value as we approach the point along any path lying entirely 
within the region of continuity. 

The definition of summability used in the discussion of the triple series 
is analogous to that of Cesaro for the simple series. Designating by simn 
the sum of the 7mn terms of the triple series > aim», lying in a rectangular 
parallelepiped 7 terms high, m terms broad, and m terms deep, in the upper, 
left-hand, forward corner of this series, and forming 


P(r +1—i) M(r+m—j) 


r(i+r) r'(m+r) I'(n+r) 


6 (r) 
(2) ym 


* Presented to the Society, December 28, 1923. 
jAnnals of Mathematics, ser. 2, vol. 24 (1922), pp. 141-166. 
{A discussion of this case has been given by Professor ©. N. Moore (cf. Mathe- 
matische Annalen, vol. 74 (1913), pp. 557-572. 
313 


314 B. M. EVERSULL 


we say that the triple series is summable (C7) to the limit 


(7) 
(3) lim Simn_ 


Imn 
and has the value of that limit, providing such a limit exists. 


I. DISCONTINUITIES ALONG PLANE SURFACES 
The triple Fourier series corresponding to the function /(., y,z) may 
by written in the form 


1 E(1/m)+ EU/n) 73 


. | | St 2’) Pimn (x, y, 2, 2", y', 2) ax’ dy’ dz’, 


where 
Pima (x, y, 2, 2", y', 2’) 
= cos [((J—1)(2’ cos [(m—1)(y’—y)] cos [(n —1)(2’—z)] 


and E(w) represents the largest integer contained in w. 

We wish to investigate first the behavior of the series at a point of 
discontinuity such that all other points of discontinuity in its neighborhood 
lie on a plane passing through that point and such that /(., y, 2) approaches 
a definite value as we approach the point of discontinuity from either side 
of the plane. Several lemmas are necessary before we can prove the 
fundamental theorem. The first two are merely stated, their proofs appearing 
in a previous paper by the writer.” 

LemMMA 1. Let BR be a region in space lying within the cube whose sides 
are + (a—o,), 8B = t(a—G), y = +(4—O01), and outside the 
sphere of radius @2 whose center is at the origin, where eo, and @2 are two 
positive constants whose sum is less than a. Then if (a, Bx) is finite 
and integrablet in R, the limit 


sin-/@ sin" 27 

lim {J y(a, B, 7) de dg dy 
lmn JR 


1, m,n —> 00 sin?« sin®8 sin®y 


will exist and be equal to zero. 


* Loc. cit., pp. 154-156, 157-158. 
+ Here and elsewhere thoughout the paper, when it is said that a function is integrable, 
it is meant that it has an integral according to the definition of Lebesgue. 
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LemMMA 2, Let R be a region within the cube described in Lemma 1, and 
ach that the point « = B = y = 0 lies within or on the boundary of R. 
Then if p(@,B,7) ts finite and integrable in R, and if 


lim g(e«,8,7) = 0, 


the limit 
1 sin? ny 
JR 


Lmn>o (Mn, sin?« sin?f sin®y 


will exist and be equal to zero. 
Lemma 3. If g,h and k are positive numbers less than n, the limit 


*h *k 9 9 
1 sin?mA sin?ny 
in —— dadBdy 
lmn 


whe 
l,m,n—> 0 o Jo Jo lin’e sin’? sin*y 


andl the seven others in which the limits of integration are respectively 
—qgtc0, Otoh, Otok; Otog, —hto0, Otok; O tog, O toh, 
—k to 0; —g to0, —h to0, 0 tok; —g to 0, 0 toh, —k to 0; 
0 to g, —h to 0, —k to 0; —g to 0, —h to 0, —k to 0, will each exist 
and be equal to {. 

The proof for the first limit follows easily from Lemma 1 and a result 
due to Fejér.* The succeeding seven cases may be reduced to the first 
by a suitable change of variable. An analogous proof is given under 
Lemma 5 of the writer’s previous paper. 

LemMA 4, If g,h and k are positive numbers less than m and we divide 
the parallelepiped whose sides are a = +g, B = th, y = tk into two 
parts R, and Rs by passing a plane through the diagonal joining two 
opposite vertices, the limits 


1 sin?Za@  sin®?ny 


will exist and will each equal >. 

Suppose the plane passes through the diagonal joining the vertices (g, h, k) 
and (—g,—h,—k). In general this plane will cut through two pairs 
of faces. Suppose for the sake of definiteness that we assume it cuts the 


(6) 


*“Mathematische Annalen, 58 (1904), p. 55. 
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faces « = +g,8 = +h. We shall chose A, as the region included 
between the dividing plane and the planes «>= +g, B= th, y= —k, 
and show that the first limit in (6) exists and is equal to >. It will then 
follow that the second of these limits exists and is equal to =, as we 
know that the sum of the two expressions in brackets approaches the limit 
unity as 7, m and n become infinite. 

The equation of the plane may be written in the form y = pa+q8 
where p and q are constants. The quantity in brackets in the first of 
expressions (6) may then be written 


1 poh 70 2g wh 2g 20 

lm n tf, | k e/J0 e/J0 eJ/0 J0 —k 
20 20 apea+-gB 20 ah 

—gvJ0 —k 


where the quantity under the integral sign in each case is the same as 
that under the integral sign in the expression (6). 

The first term of (7) approaches + as a limit as 7, m and n become 
infinite, by Lemma 3. If we make the substitution « —= — «’, B= — 8, 
y =-—-/ in the third and fourth terms, the third combines with the fifth 
to form an expression which by Lemma 3 approaches + as a limit as 
1, m and m become infinite, and the fourth combines with the second to 
form an expression which approaches 4 as 7, m and n become infinite, by 
the same lemma. Hence in this case the limit of the first of expressions (6) 
exists and is equal to }. 

The plane through the diagonal joining the vertices (g, h,k) and 
(—g,—h,—k) may also be perpendicular to one of the planes, say 
to y =k. In this case we take R, to be the region included between the 
plane through the diagonal and the planes « = g, B=—h, y= tk. 
The quantity in brackets in the first of expressions (6) will then be 


1 (70 20 thalg ek 


where the expression under the integral signs is the same as in (6). From 
1 


Lemma 3 the first term of this expression approaches ; as 7, m and n 
become infinite. If in the third term we make the change of variable 


axe’, 
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this term combines with the second to produce an expression which, from 
Lemma 3, also approaches } as 7, m and n become infinite, and hence for 
this case also the first limit in (6) exists and is equal to >. 

Hence this is true for all positions of the plane through the diagonal, 
and as shown before, the second limit in (6) must also exist and be equal 
to : for all positions of the plane. The lemma is therefore proved. 

Lemma 5. If R is the region described in Lemma 2, the limit 


1 sin?Za@ sin? sin? ny 
lim | | dadBdy 
lmnJ J Jr sin*y 


1 le sin? sin® ny 
J JR 


will exist and be equal to zero. 
The expression in brackets may be written in the form 


ine whe 

sin*J/@ sin? ny 

« sin*8 sin* y 
It may easily be seen that the product of the three terms in parentheses 
under the integral sign satisfies the condition imposed on ¢(a@, 8, y) in Lemma 2, 
and hence the lemma follows at once. 

Lemma 6, Let R be a region within the cube described in Lemma 1, and 
without the cube where is an 
arbitrarily small positive quantity and y(«@, 8,7) a function that is integrable 
mk and remains finite in that part of R which includes all points whose 
codrdinates satisfy one or two of the conditions 


(8) —ecace, 
Then the limit 


(9) lim de dB dy 


1,m,n—>o e 


) sin?2a@ sin?m~f n 
e e 


will exist and be equal to zero. 
1e quantity in (9) may be written in the form 
The quantity in (9) may b tt the f 


40 SS Set SS 


— 
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where the integrands are the same as that in (9), A’ that part of R which 
includes all points whose coérdinates satisfy one or two of the conditions (8), 
and R” the remainder of R. 

By Lemma 1 the first term of (10) approaches zero as a limit as 1, m 
and » become infinite. 

Consider now the second term of (10): 


Y 
R’ ( 7 ) sin? « sin? 8 sin? y B | 


| 
Imnsin® g(a, B, y)| de dB d; 


where 7 is the smaller of the quantities « and d, 6 being the least distance 
between and the cube ta, B= +n, y= +a. The right-hand 
side of this inequality may be made as small as we please by choosing /, m 
and » sufficiently large and therefore the second term of (10) approaches 
zero as /,m and » become infinite. Hence the limit (9) exists and is equal 
to zero and the lemma is proved. 

Before proving the next theorem we shall define an expression which 
we wish to use in the statement of this and succeeding theorems. A critical 
region associated with the point (x,, y¥,,%) is defined as a region made up 
of points whose coérdinates satisfy one or more of the inequalities 


1 fff sin?7@ sin?mf sin®?n 
e e 


lmn,. 


where ¢ is an arbitrarily small positive quantity. 

We are now ready to prove the following theorem: 

THEOREM I. Jf f(x, y, 2) is integrable in the region 


(11) 


and (2,, 41, 2) is a point of discontinuity of f(x, y, 2) such that every 
other point of discontinuity in the neighborhood of (2,, y:, 2) lies on a plane 
passing through that point, and the function approaches a definite value as 
we approach the point from either side of the plane, the series (4) will be 
summable (C1) at that point to a value half way between the limiting values 
of the function, provided f(x, y, 2) remains finite in some critical region 
associated with the point (1%, y1, %). 
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For the series (4) at the point (2,, ¥,, 2,) we have 


1) 7, 
Simn(%15 _ | 
lImn Simnz*® J 


") 


where Sinn is defined by equation (1). Making the change of variable 


(13) a’ 2 «a, y'— — 28, 
we have 


S(mt2e@, 


» 


(14) 
“Ine , 
sin*f sin*y 

In order to prove the theorem we must show that as 7, m and x become 
infinite, the right-hand side of the expression (14) is equal to one-half the 
sum of the limiting values of the function. 

Since all the points of discontinuity in the neighborhood of (2,, ¥,, 2; ) lie in 
a plane, we can determine a rectangular parallelepiped of dimensions 24, 
2h, 2k whose sides are parallel to the coérdinate planes and such that 
the point (a, ¥%:, 2,) lies at the center of the parallelepiped and all other 
points of discontinuity lie on a plane through (2, y;, 2,) either parallel to 
two opposite faces of the parallelepiped or passing through one of its 
diagonals. 

If we make the transformation (13) this becomes a rectangular parallelepiped 
whose center lies at the point «e = 8 = y = 0. We may call this region R’, 
designating the remainder of the region of integration in (14) by R”. We 
then have 


y(1) > > > 
lmn J J Je J Jr’ 


where the quantity under the integral signs is understood to be the same 
as that in equation (14). 
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By virtue of Lemma 6 the second term of the right-hand side of this 
equation approaches zero as J, m and n become infinite. 

The region R’ is divided into two parts, Rj and R3, by the plane of 
discontinuity and hence the first term on the right-hand side of (15) may 
be written in the form 


(16) lmnn® | | J Jr,’ 


the expression under the triple integral signs being the same as in 
equation (14). 

If we represent by f, and f2 the values which /(, y, z) approaches as 
we approach (2, ¥%,, 2,) through the regions Rj and FR; respectively, we 
may write the first term of (16) in the form 


1 sin?mf sin?ny 
J g(«, 8,7) sin?8 sin? y dedBdy 


sin*y dedBdy, 


sin*/« sin?mf sin’ ny 
J JR; 


where g(«, 8, vy) = f(m It then follows 
from Lemma 2 that the first term of this expression approaches zero as 
1, m and m become infinite. The second approaches 3f, as 7, m and n 
become infinite, by virtue of Lemma 3 if the plane of discontinuity is 
parallel to two opposite faces of the parallelepiped, or of Lemma 4 if it 
passes through a diagonal. 

Similarly it may be shown that the second term of (16) approaches if, 
as /, m and nm become infinite and therefore (16), and hence the left-hand 
side of (14), approaches 3(f, + /2) as 1, m and » become infinite and the 
theorem is proved. 

For points on the boundaries of the region (11) except at the edgés and 
vertices, the results are analogous to those for points in the interior and 
may be obtained by methods similar to those of Theorem I. These results 
are given in the following corollary: 

CoroLLaRy. Jf 2, and y, are values of x and y lying in the interval 
(—a and if f (2, y, is integrable in the region 
(11) and furthermore is such that the limits 


(17) lim y; 2), lim z) 
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evist, the Fourier development of f(x, y, 2) will be summable (C1) at the 
points (215 and (x1, —7) and its value will be one-half the sum 
of the two limits (17), provided f(x, y, z) remains finite in some critical 
regions associated with (7,41, and 41, —m). Analogous statements 
with corresponding conditions may be made for points on the other faces. 

We shall next consider some of the results obtained when (7, y,, 2,) is 
a point of discontinuity of f(x, y, z) such that all other points of discon- 
tinuity of the function in the neighborhood of that point lie on two planes 
through the point. When the planes are parallel to two of the coérdinate 
planes the results are simple; this case, which is of special importance since 
itis analogous to the case where the points lie on the edges of the region 
of periodicity, is considered in the following theorem: 

THEOREM I. Jf f(a, y, 2) is integrable in the region (11) and (a, m1, 2%) 
is a point of discontinuity of f(x, y, 2) such that every other point of dis- 
continuity in its neighborhood lies on one of two planes which are parallel 
to two codrdinate planes and pass through (a1, 1,2), and if furthermore 
the function approaches a definite value as we approach the point through 
each of the four regions into which the planes of discontinuity divide the 
neighborhood of that point, the series (4) will be summable (C1) at (x1, %1, 21) 
to a value which is one-fourth the sum of the four limiting values of the 
function, provided f(x, y, 2) remains finite in some critical region associated 
with the point (a, 21) 

The proof of this theorem is similar to that of Theorem I. The follow- 
ing corollary may also be established by similar methods: 

CoroLLARY. Jf a, is a value of x lying in the interval (—na <x <2) 
and if f(x,y, 2) is integrable in the region (11) and furthermore is such 
that the four limits 
(18) lim S(x,y; 2) 


exist, the Fourier development of f(x, y, 2) will be summable (C1) at each 
of the four points (a,, +m, +m) and to a value which is one-fourth the 
sum of the limits (18). Similar statements may be made under analogous 
conditions for points on the other edges. 

temark.—The conclusions of each of the first six lemmas hold equally 
well if instead of becoming infinite together, 7, m and m become infinite 
in any other manner. Hence Theorems I and II and their accompanying 
corollaries hold however J, m and m become infinite. 

In all other cases where the point of discontinuity is such that all the 
points of discontinuity in its neighborhood lie on two planes through the 
point, the value which the expression 
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(19) lmn 


approaches as 7, m and » become infinite depends on the positions of the 
planes and also the manner in which 7, m and m become infinite. Before 
considering some of these cases it will be necessary to obtain some 
preliminary results. 

For the purposes of our work we may divide the parallelepiped of 
Lemma 4 into four parts in two ways, namely by two planes each of which 
passes through a diagonal joining two opposite vertices of the parallelepiped 
or by two planes one of which passes through a diagonal and the other 
passes through the center and is parallel to one of the faces. We may 
now state the lemma: 

Lemma 7. (a) If we divide the parallelepiped of Lemma 4 into four 
parts R,, Rz, R;, Ry by either of the above methods, and if each of the 
ratios m/n, n/l, l/m approaches a finite limit or becomes infinite as 1, m and n 
become infinite, the limits 


1 sin2le sin? mB sin?ny 
(202) lim —— — daddy (i = 1,2,3,4) 
l,m,n—> 0 lmn J J. Rk, Sin 8 


will exist and will each approach a limit between 0 and a*/2. Moreover 
the sum of two limits in which the regions of integration are adjacent is 


m*/2, and the limits for which the regions of integration are opposite are 
equal. If any of the ratios between 1, m and n oscillates, the limits will 
oscillate between 0 and 2*/2, their values being subject to the same restriction 
as above. 

(b) If the intersection of the planes is parallel to the «- (B- or y-) axis 
the results are the same as above except that only the ratio m/n (n/1 or l/m) 
affects the values of the limits, which are fixed if this ratio remains fixed 
and oscillate if it oscillates. 

We shall consider, first, part (a). That the value of each limit lies 
between 0 and 7°*/2 or is equal to one of them, subject to the restriction 
that the sum of two limits in which the regions of integration are adjacent 
is equal to 2*/2, is obvious from Lemma 3 if the plane used as the dividing 
plane is parallel to a coérdinate axis, or from Lemma 4 if it is not. 

That the limits are the same for opposite regions of integration may 
be seen by making suitable changes of variable of the type « = —e’. 

Each of the limits (20) depends on a limit of the type 


1 sin?’mf sin? ny da dp dy 

L.mn—>o lmn Jo Jo Jo sin?« 
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which, if we apply Lemma 5 and make the change of variable 


la=«a', mB=f, ny =y', 
reduces to the form 


‘ sin*e sin*y 
This limit depends on the values of the ratios between /, m und m and 
hence when any of these ratios oscillates the value of the limit oscillates. 
To prove the lemma for case (6) we may assume that the intersection 
of the planes is parallel to the «-axis. Each of the limits (20) depends on 
a limit of the type 
ah . . 
1 sintZa sin?’mf sin? ny 
lim 2 da dp dy 
imn>olmnJo Jo Jo sin°f sin*y 


which reduces to the form 


(99) sin 2a ré_sin’y dy |\dBi de 


by Lemma 5 and the change of variable (21) if we set 


r= lim —k. 
mn—>o 
The proof then follows as in part (a). 

In this case we may obtain an exact value for the limits (20) since 
they depend on triple integrals of the type of (22), which we may evaluate 
without great difficulty. The double integral in brackets in (22) has been 
evaluated by C. N. Moore*, and applying his results we have for the value 
of (22) 

ne 4 { (r +1) log (r+1)— 2r logr + (r—1) log (r—1) 

8 2 | 4 
(r—1) log(r—1)—(r+1) log(r+1) 1 
4r r 


1 1 
+. + +-lo>n, 


16 (r = 1), 
2 4 


log (irr) (r-+ 


1 


5 


(0<r<1). 


* Loc. cit. p. 564. 
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We may now consider the summability of f(z, y,z) at a point of dis- 
continuity such that all other points of discontinuity in its neighborhood 
lie on two planes both of which are not parallel to a coérdinate plane. 
The proofs may easily be carried out by methods similar to those of 
Theorem I except that Lemma 7 is used instead of Lemmas 3 and 4. 

Let f(2, y,z) be a function that is integrable throughout the region (11) 
and (2, #1, 2;) a point of discontinuity of f(a, y, z) such that all the dis- 
continuities in its neighborhood lie on two planes which pass through that 
point and whose intersection is not parallel to a coérdinate axis, and further- 
more let f(a, y, 2) remain finite in some critical region associated with 
(23, %1, 2) and approach a definite value as we approach (2, y, z) through 
each of the regions into which the planes of discontinuity divide the neighbor- 
hood of that point. Then if 7, m and m become infinite in such a manner 
that any of the ratios n/m, l/n, m// oscillates, the value of (19) will oscillate 
between and where fi, fs, fa are the values which 
the function approaches, assigned in rotation. If the ratios become infinite 
or remain finite the expression (19) will approach a finite limit between 
‘(fit fs) and i(fo+f,) or equal to one of them. 

If f(x, y, z) and (7,, %, %) satisfy the conditions of the preceding para- 
graph except that the intersection of the planes of discontinuity is parallel 
to one of the coérdinate axes, the only ratio which affects the limit of the 
expression (19) is n/m, l/n, m/l, according as the intersection of the planes 
is parallel to the a-, y- or z-axis. If this ratio oscillates the value of the 
limit (19) will oscillate between i(f,+.f,) and $(f.+-/1), and if it ap- 
proaches a finite limit or becomes infinite (19) will approach a limit between 
these values or equal to one of them. The exact value may be obtained 
by determining the values of the limits (20) from the expressions for the 
integral (22). 

In particular we may derive the following theorem: 

THEOREM III. Let f(x,y, 2) be a function which is integrable in the 
region (11) and (a, y%:, 2,) a point of discontinuity of f(x, y, 2) such that 
every other point of discontinuity in its neighborhood lies on two mutually 
perpendicular planes which pass through (2,, 41, 2) and whose intersection 
is parallel to one of the codrdinate axes. Then if f(x, y, z) remains finite 
in some critical region associated with (a1, y:, 2) and approaches a definite 
value as we approach that point through each of the four regions into which 
the planes of discontinuity divide its neighborhood, the series (4) will be 
summable (C1) at (a1, %1, 2) to one-fourth the sum of the four limiting 
values of the function, provided the ratio n/m, W/n, or m/l, according as the 
intersection of the planes of discontinuity is parallel to the x-, y-, or z-axis, 
approaches unity as l, m and n become infinite. 
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Results similar to those obtained in the above discussion hold where the 
discontinuities are distributed on a broken plane instead of on two inter- 
secting planes. 

We may also carry through discussions of the cases where the dis- 
continuities lie on three or more planes intersecting at a point of dis- 
continuity, but as these results are in general complicated we shall confine 
ourselves to stating one particular case from which we may deduce the 
behavior of the series at the vertices of the region (11). The results are 
contained in the following theorem and its corollary, the proofs being carried 
out by a method similar to that used in Theorem I: 

THEOREM IV. Jf f(x, y, 2) is integrable in the region (11) and (2, #1, 21) 
is a point of discontinuity of f(x, y,z) such that every other point of 
discontinuity in its neighborhood lies on three planes through that point and 
parallel respectively to the coirdinate planes, and if moreover f(x, y, 2) 
approaches a definite value as we approach (2, 4, 2) through each of the 
eight regions into which the planes divide the neighborhood of that point, 
the series (4) will be summable (C1) at the point (2,1, #1, 2) and to a value 
which is one-eighth the sum of the eight limiting values of the function, provided 
flix, y,2) remains finite in some critical region associated with the point 
(11,415 

CorouLaRY. If f(x, y,z) is integrable in the region (11) and if the eight limits 
(23) lim (ayy, 2) 

erist, the series (4) will be summable (C1) at the eight vertices of the region 
(11) and its value will be one-eighth the sum of the eight limits (23), provided 
the function remains finite in some critical region of each of the eight vertices. 


IJ. DISCONTINUITIES WHICH LIE ALONG CURVED SURFACES 

We shall now consider the behavior of the series at a point of discontinuity 
such that all other points of discontinuity in its neighborhood lie on a curved 
surface through the point and such that the function /(2, ¥, 2) approaches 
a definite value as we approach the point of discontinuity from either side 
of the surface. This surface is assumed to have a tangent plane at the 
point of discontinuity under consideration and to be of such nature that 
any plane through the point intersects the surface in a finite number of 
curves or not at all. 

Before proving our main theorem we shall have to prove several lemmas. 

LemMMA 8. The integral 


sin? a sin®f sin’ 
(24) | | | — dadBdy 
eJ/0 2/0 a 


326 B. M. EVERSULL (July 


converges for all values of 4 and w and represents a function of 2 and p, 
w(d, w), which is continuous for all values of 4 and » and approaches n*/8 as 
a limit as 4 and w become infinite in the same or opposite directions. 

We shall consider first the case where 4 and mw are both positive. If 


we let 
va pp 


sin? y 
pw, a, 8) = —a dy 


we may write the integral (24) in the form 


(25) @, B) dad. 
Jo J0 a” 


It may readily be seen that (A, «, 8) is a continuous function of 2, «, 8 
for all values of these arguments and that it remains finite for all such 
values. Therefore the integrand of (25) is a continuous function and as the 
double integral (25) is uniformly convergent, w(4, ), the function repre- 
sented by (25), defines a function of 4 and w which is continuous for all 
values of 4 and p. 

If A and w are of opposite sign, say 4 is negative, we may make the 
substitution « = —«’ and the proof follows in the same manner as before. 
If both 2 and mw are negative, we set « = — a’, 8 = — #’, and the proof 
follows similarly. 

LemMA 9. The integral 


sin?a sin? sin®y 


ce B* 


(26) -dadBdy, 


where w remains fixed, converges for all values of 4 and represents a function 
of 2, 2(4), which is continuous for all values of 4 and remains finite as 4 
becomes positively or negatively infinite. 
Defining 
sinty 
df dy 


by z(4, «), we may write the integral (26) in the form 


ws sin? « 
(27) | x(A4, da, 
e/ ( 


) 


A 
——« 
J0 J0 
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It is obvious that y(4, «) is a continuous function of 4 and « for all values 
of these arguments and that it remains finite for all such values. Hence 
the integrand of (27) is a continuous function of 2 and « for all values 
of these arguments and since the integral is uniformly convergent for all 
values of 4, it represents a function, 2(4), which is continuous for all 
values of 4. 

As 24 becomes infinite 4(4, «) approaches the limit 2*/4, provided « + 0, 
and hence 2(4) approaches 2*/8 as a limit and therefore remains finite for 
all positive values of 4. 

For the case where 4 is negative, the proof may be carried out by 
making the change of variable « = —@’, 

LemMA 10. If y=G(e@,8) represents a curved surface whose tangent 
plane at the origin is y =4a-+ mB and whose intersection with the «,B 
plane is B=g(ea), and if y=—G(a,8) is intersected in only a finite 
number of curves in the neighborhood of the origin by any plane passing 
through that point, the limits 


ah pk . 
lim J J P sin? sin’ 8 sin? y dadpa 
imn a ¢ dy 
l,m,n—>o lmn 0 J0 0 a B 


(28) 
vk 
¢ sin* 28 sin® 
lmn B 


g(@) Ga, in? « sin? sin® 
lim 2 i J, 


aA + 


where h and k are constants between 0 and mn, will exist and be equal to zero 
provided that 4 and w are not zero in (28) and that w is not zero in (29).* 
We shall carry out the proof where 4 and w are both positive, the other 
cases following by a change of variable as in Lemmas 8 and 9. 
Consider first the expression (28), writing the quantity in brackets in 
the form 


*If the equation of the surface or tangent plane is such that it can not be put in the 
form used in the statement of the lemma, we may use « or # as the left-hand side of the 
equations, modifying the remainder of the statement and the proof accordingly. 
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Imn l. 0 J0 Jo 
Imn \Jo Jo Jo 0 J0 Jo J 
ak phat 


where 0 < p< h,0<q<k, the expression under the triple integral signs 
being the same as that in (28). It is obvious that the second and third 
terms of (30) approach zero as 7, m and m become infinite. Then if we 
can show that p and q may be so chosen that for sufficiently large values 
of 7, m and » the first term of (30) may be made as small as we please 
in absolute value, the lemma will have been demonstrated for the 
expression (28). 

Assume for the sake of definiteness that the portion of the surface above 
the «, 8 plane lies above the tangent plane until possibly it intersects it. 
We may choose a 4’ and w’ greater respectively than 4 and mw, such that 
= and such that the plane y = intersects the surface 
in a curve whose projection on the «, A plane lies entirely within the 
projection of any curve in which the tangent may cut the surface in the 
octant O—YYZ. If we draw a rectangle two of whose sides lie along 
the « and 8 axes, within the area included between these axes and the 
projection of the intersecti and y = @(a, 8) in this 
octant, the vertex opposite the origin being (p, qg,0), we have 


lImnJo Jo Jo lmn Jo Jo Jo lmn Jo Jo Jo 


the expression under the integral signs being the same as that in (28). 
Making the transformation (21) we see that, p and q being fixed, we can 
choose / and m so large that the first and third members of the inequality (31) 
become as near and w(nd/l, nu/m) respectively 
as we wish, this function being defined as in Lemma 8. But since by that 
lemma this function is continuous, it follows that for 2’ and w’ sufficiently 
near to 4 and w respectively, and w(nd/l, nu/m) 
differ for all values of 7, m and » by as small a quantity as we please.* 


* The restriction 4+0, # +0 is of course essential for this conclusion. When these 
conditions do not hold, we must limit the freedom of variation of 1, m and as is pointed 
out in Lemma 11. 


[July 
328 


1924] TRIPLE FOURIER SERIES 329 


Hence from (31) p and g may be so chosen that for 7, m and » large 
enough the first term of (30) is as small in absolute value as we please, 
and the lemma is proved for the expression (28). 

To prove it for (29) we write the quantity in brackets in that expression 
in the form 


+ 


l m n Jo 0 2/0 eJ0 e/0 
oh og (a) eG (a, 8) og (a) -G(a, 8) 
lmn |Jo Jo 0 J0 J0 


A 


sh 


what pp > ft whe 4 
+7. n J, j P 


where 0 < +r <h, the expression under the triple integral signs being the 
same as that in (29). 

Here as in (30) the second and third terms approach zero as 7, m and n 
become infinite. Then we have only to show that 7 may be so chosen 
that for sufficiently large values of 7, m and » the first term of (32) may 
be made as small in absolute value as we please, and the lemma will have 
been proved for the expression (29). 

We shall again use the surface which we used in the proof for the 
expression (28), considering that part which lies in the octant 0—X 1’ Z. 
We may then choose a 2’ greater than 4 and such that the plane y = 2’« +8 
intersects the surface in a curve whose projection on the «, # plane in 
the octant O—X1"Z lies entirely within the projection of any curve in 
which the tangent plane may cut the surface in that octant. Then if r is 
the « coérdinate of the first point to the right of the origin in which the 
tangent plane intersects the surface, we have 


a- G(«, 
+ pp 
Imn 


where the expression under the integral signs is the same as before. 


— A 
| 
(32) + | 
(33) 
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Making the transformation (21) we proceed by a method analogous to 
that used in the proof for the expression (28), using Lemma 9 instead of 
Lemma 8, and the lemma is proved for the expression (29). 

LemMA 11. Jf all the conditions of Lemma 10 are satisfied by y = G(«, 2) 
except that the tangent to the surface passes through the «- or B-axis or 
coincides with the «a, 8 plane we have the following results: 

(a) if the tangent plane passes through the a-axis (2 = 0), the limit (28) 
exists and is equal to zero provided n and 1 become infinite in such a manner 
that 
(34) nil < K, 


where K is a positive constant. Similarly if the tangent passes through 
the B-axis (uw = 0), the limit (28) exists and is equal to zero provided 


that n and m become infinite in such a manner that 


nim K; 


(b) if the tangent plane passes through the a-axis the limit 


1  sinte sin?B sin*y 
lmn Jo g(@) J 0 a B 7 


will exist and be equal to zero provided l, m and n become infinite in such 
a way that the ratio between any two of them remains finite; 
(c) if the tangent plane coincides with the a, B plane the limit 


/ k . 9 
lim 1 | sin?@ sin®y le dB d 
wo lmn Jo Jo Jo 


l,m,n—> 


will be exist and be equal to zero provided 1, m and n become infinite as 
provided in (b). 

The proofs involved in this lemma may be obtained by making slight 
changes in Lemma 10. For example, in the first part of (a) where 4 = 0, 
we cannot say that if 2’ is taken near enough to 4, w(mnA’/1, np’/m) 
and w(n4/l, n/m) differ by as small a quantity as we please for all 
values of 7, m and n. The statement is true for only such values of / 
and n as satisfy the inequality (34). The remaining proofs may be carried 
out by making similar changes in Lemma 10. 

Before proceeding with the next theorem we shall define a term which 
we wish to use in the statement of the theorem. If the quotient in (3), 
corresponding to a triple series > aimn, approaches a limit as 7, m and 
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become infinite provided the ratios between /, m and m remain less in 
absolute value than a positive constant, we say that the triple series is 
vestrictedly summable (Cr) and that its limit is the value approached. If 
only one ratio needs to be restricted we say that the series is restrictedly 
summable with respect to the quantities involved in that ratio. 

THEOREM V. (a) Let f(x,y,z2z) be a function which is integrable 
throughout the region (11) and (1,1, 21) a point of discontinuity of 
flv, y, 2) such that every other point of discontinuity in its neighborhood 
lies on a surface satisfying the following conditions: (1) the surface has 
a tangent plane at the point (a1, y1, 21) which ts not parallel to any of the 
coirdinate axes; (2) no plane through the point (a1, y1, 21) intersects the 
surface in an infinite number of curves in the neighborhood of that point; 
and (3) the function f(x, y,2) approaches a definite value as we approach 
(11,41, 21) from either side of the surface. Then the series (4) will be 
summable (C1) at the point (a1, y1, 2) to a value half way between the 
limiting values of the function at that point, provided f(x, y, 2) remains 
finite in some critical region associated with (a1, y1, 21). 

(b) If all the other conditions are fulfilled, but the tangent plane to 
the surface at (x1, 41, 21) is parallel to a codrdinate axis or plane, the 
series will be restrictedly summable to one-half the sum of the limiting values 
of the function at (a1, y1, 21). If, however, the plane tangent to the surface 
at (x,y,z) is parallel to only one of the codrdinate axes and remains 
tangent to the surface in a line parallel to that axis, only one ratio needs 
to be restricted, the ratio being nim, l/n, m/l, according as the line of 
langency is parallel to the x-, y-, or 2-axis. 

We shall consider first the case (a). We may choose a rectangular 
parallelepiped R with (2, y:, 2) as its center, such that the tangent plane 
passes through one diagonal and all points of discontinuity in the parallele- 
piped lie on the surface of discontinuity, and small enough so that the sur- 
face of discontinuity does not intersect the plane tangent to the surface within 
the parallelepiped, except perhaps at the point of tangency. 

The Fourier series corresponding to f(x, y, z) may then be written 


(1) 
Stmn (x7, a) 1 
lmn lmnn*. i+ 1+27) 
sme 
sin*la sin*n 
da dp dy 
sin*« sin? y 


| t2e,n +22, a+2y7) 


lmnn* 


sin?2a@ sin?mA sin?ny 
sin* « sin® 8 sin*y 


24 
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where R’ represents the rectangular parallelepiped into which F is trans- 
formed by the change of variable (13) and R” is the remainder of the region 
of integration of the integral (14). The second term of (35) vanishes as 
l, m and n become infinite, by Lemma 1. Then in order to prove our 
theorem we have only to show that as J, m and m become infinite the 
first term approaches a value which is one-half the sum of the limiting values 
of f(x, y, z) at 

Let the surface of discontinuity divide the parallelepiped into two parts, 
Ri and R; and let /f; and /f2 be the limiting values of the function as we 
approach the point of discontinuity through A; and #3 respectively. Then 
the first term of (35) may be written 


where the quantity under the triple integral signs is the same as in (35), 
We shall show that the first term of (36) approaches ;/, as a limit as /, m 
and » become infinite. 

The tangent plane to the surface at (z,, y,, 2,) divides the region Rj into 
two parts; let D; be that one which lies on the same side of the tangent 
plane as Rj does of the surface of discontinuity. Then by Lemma 10 


1 sin?’mf sin?n 


sin*la sin® in? 


approaches zero as 7, m and m become infinite, and hence by Lemma 5 


sin* sin sin? ny 


J, sin*/« sin?m8 sin® 
rR, sinta sin?B 


approaches zero as 7, m and n become infinite. But by Lemma 4 the first term 
of (37) approaches 2*/2 and hence the second term must also approach 2°*/2 
as 1, m and m become infinite. Therefore the first term of (36) approaches ;/; 
as 1, m and n become infinite. Similarly it may be shown that the second 


(37) 
(38) 
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term of (36) approaches the value 3/; and therefore the whole expression (36) 
and hence the right-hand side of (35) approaches as 7, m and 
become infinite, and part (a) of the theorem is proved. Part (b) may be 
proved similarly by using Lemma 10 instead of Lemma 9. 

We wish now to show that restricted summability is all that may be 
established for part (b). Since the proofs are analogous for the cases 
where the tangent plane is parallel to a coérdinate axis or to a codrdinate 
plane, we shall give the proof only for the latter, considering the develop- 
ment of a function that is continuous throughout the region (11) except 
on a curved surface lying above the x, y plane and tangent to it at the 
origin, and that approaches f, as we approach the surface of discontinuity 
from above the surface and /2 as we approach it from below. 

From Moore’s generalization of Fejér’s theorem, 


00, 00 
((0',0,0)= f(z’, y', 2’) 


-cos[(m—1)y'] cos[(n —1)2’] dx’dy'dz’ +0). 


The series on the right may be integrated term by term if we multiply it 
by a continuous function of x, and therefore 


nat 
7 fs (x 0, 0) cos[(/ ] dx 


m=1,n 


cos[(m —1) y’] cos[(n —1)2’] dz'dy'dz’, 


and by the generalization of Fejér’s theorem, 


|_| 


- cos[(7 —1)z’] cos[(m —1)y’] cos[(n — 1) 2’] ae 


24° 
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and the expression (1/7 m2) [ Sion (0,0,0)] approaches /, if we let m and » 


become infinite holding 7 fixed, and then let / become infinite. If we let 
1 and n become infinite holding m fixed and then let m become infinite the 
series is again summable to /2 at the origin, but if we let / and m become 
infinite and then mn, the value to which the series is summable at the 
origin is +(f,+ 2) and hence the triple series cannot be summable in 
the ordinary sense. As pointed out above, it may be shown in analogous 
fashion that restricted summability is all that can be proved when the 
tangent plane is parallel to one codrdinate axis. 

By methods similar to those of Theorem V and its preliminary lemmas 
we may obtain the value to which the function is summable at points of 
discontinuity such that all other points of discontinuity in their neighborhood 
are on two or more curved surfaces through the point, or lie at the inter- 
section of two or more such surfaces. The results in these cases are 
analogous to those of Theorems II, IJI and IV, except that when the 
tangent to one or more of the curved surfaces is parailel to a codérdinate 
axis we have restricted instead of ordinary summability. 

UNIVERSITY OF CINCINNATI, 

CrxcrnnaTI, OHIO. 


AN UNUSUAL TYPE OF EXPANSION PROBLEM* 
BY 


M. H. STONE 


In the majority of papers on the expansion problems associated with 
linear differential systems or with integral equations, the chief aim is a 
proof that a function, arbitrary within the limits of certain restrictions 
as light as possible, can be expanded in a uniformly convergent series of 
functions arising from the differential system or the integral equation. As 
exceptions to the rule we may cite a number of papers dealing with dif- 
ferential systems for which the boundary conditions are of irregular type.? 
In the following paragraphs we present a study of the linear differential 
system of the first order, without the customary limitation that a certain 
coefficient in the differential equation remain positive; the results obtained 
are in sharp contrast not only with those found for the usual expansion 
problems and those already mentioned in connection with irregular boundary 
conditions, but also with the facts for the precisely analogous second order 
differential system.t A second point of some interest is the appearance of 
a normal orthogonal set of trigonometric functions differing from the nor- 
malized Fourier set but as closely related to it as are the sine and cosine 
sets. It is interesting that this exceedingly simple differential system seems 
to have escaped attention entirely. So far as we know, the existence of 
this gap in the theory of boundary value and expansion problems was first 
pointed out by Professor Birkhoff in a recent course of lectures on the 
theory of linear differential equations. 


I. THE FORMAL EXPANSION PROBLEM 


We shall study here the two linear differential systems 
(1) u’ +Apu = <2z<1; u(0)—u(1) = 0; 


(2) u’ u(0)+u(1) =0; 


* Presented to the Society, March 1, 1924. 
j Jackson, Proceedings of the American Academy of Arts and Sciences, vol.51 
(1915-16), pp. 383-417; Hopkins, these Transactions, vol. 20 (1919), pp. 245-259. 
; Mason, these Transactions, vol. 8 (1907), pp. 427-432; Lichtenstein, Rendiconti 
del Circolo Matematico di Palermo, vol. 38 (1914), pp. 113-166. 
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in which 4 is a complex parameter and p is a bounded summable function 


with Lebesgue integral [rae equal to unity, which satisfies the further 
requirement that a finite set of non-overlapping intervals @,, @2, ---, @n can 
be determined such that 
(a) the set covers the unit interval just once; 
(b) on either p> O almost everywhere (i.e., except upon a set of 
points of zero measure) or p <0 almost everywhere; 
(c) if p>O on o; almost everywhere, then p < 0 almost everywhere 
on the two adjacent g-intervals and vice versa. 
It has been customary to assume p > / > 0 on the unit interval in dis- 
cussions of expansion problems; in the present paper we dispense with 
this requirement. It will be apparent in the course of the development 
that there is no essential generalization in supposing p of the character 
already described rather than continuous with a finite number of changes 
of sign. At this point we may also remark that the more general dif- 
ferential system 


u +(Apta)u 0, <2<b, 
Au(a)+ Bu(b) = 0, 


*1 
where 4 and p are restricted as before except that I, pdzx is merely dif- 


ferent from zero, g is any summable function, and A and B are real numbers 
different from zero, can be reduced by a series of real linear transformations 
of the independent variable and of the parameter together with a real 
transformation of the independent variable in the form 


u, A(x), A(a,) + 0, 0< <1, 


to either (1) or (2) according as AB <0 or AB > 0; the transformation 
can be effected in the sense that the new differential equation is true except 
possibly on a set of zero measure. If complex transformations are ad- 
mitted (1) can be taken into (2), and conversely, in an analogous manner. 

The adjoint differential systems for (1) and (2) are readily computed as 
(3.1) Lpv —v(0)+v(1) = 0; 


(3.2) ‘tipv = 0, v(0)—v(1) = 0. 
The differential equation «’ + Apu — O has the function 


AP(@) P(x) = 
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as a solution in the sense that w satisfies the differential equation except 
possibly on a set of measure zero. That this is the case may easily be 
verified by substitution, if it is remembered that P’(2) = p(a) almost 
everywhere, in accord with the theory of Lebesgue integration.* Further- 
more, any other solution «, which is a Lebesgue integral is a constant 
multiple of «. To prove this we notice that 


| — | — 


dx u 


On integration we find that «#,/«u = a constant. In the same way, the 
adjoint differential equation —v’ + = has a solution e?, unique 
except for a constant factor. When these solutions u and v are substituted 
in the boundary conditions of the two differential systems and in the adjoint 
boundary conditions, the characteristic equations 


(4.1) 
(4,2) 


are obtained. The roots of these equations, the characteristic numbers, 
are immediately seen to be 


(5.1) Ak 2k 


(5.2) Ay = mit2kni, 


In case we had admitted the equality f pdx = Owe should have found 
here that every 2 is a characteristic value in (1), and that no 4 is a 
characteristic value in (2). Under our hypotheses the characteristic values 
form the discrete sets of numbers just determined. For solutions of (1) 
and (3.1)*corresponding to distinct characteristic numbers 4;, 4x, we have, 
from the differential equations, 


v1 
(Aj — Ax) | purjdx = + = 
e 


0 


* Lebesgue, Lecgons sur l’Intégration, Paris, 1904, pp. 124-125; de la Vallée Poussin, 
Intégrales de Lebesgue, Paris, 1916, § 67. 


| k=0, 41, 
k=0, +1, +2,---. 
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On account of the boundary conditions for (1) and (3.1) the last expression 
vanishes, so that 


= 0. 
eJ0 


On the other hand, for solutions corresponding to the same characteristic 
number we are led to 


| = ot hP@ - pdz = 
eJ0 eJ0 
In similar fashion we find for (2) and (3.2) 


J, purvjda \1 = 


By analogy with Fourier series the following formal expansions of an 
arbitrary function integrable in the sense of Lebesgue are suggested: 


+ 00 
(6.1) Ak Uk; = dx; 
oO 
+ 00 a1 
eJ0 


Uk; = pfu dx; 


J0 
+- a1 
0 


If we perform the formal operation of collecting the conjugate complex 
terms in these series, we find two related formal series of real terms with 
real coefficients. For the series (6.1) and (7.1) the terms with subscripts 
+k and —k respectively are conjugate complex; and the grouping of these 
terms yields one real series 


V2 cos2ka P(x) + By V2 sin2ka P(x) }; 


1), 

‘1 = 
A, = | pf V 2 cos2ka P(x)dz, 
0 

| pth 2sin2kaP(x)dzx. 


0 


(8.1) Ap 
1 


= 
Ay 
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Similarly, the pairs of conjugate complex terms in the series (6.2) and (7.2) 
are those with subscripts k = 0, —k—1 respectively; and the series ob- 
tained by collecting terms is 


= 
{ Ap V2 cos(2k-+1)aP(r)+ By, V 2sin(2k+1)a )} ; 


pf V 2cos(2k+1)aP(ax)dz, 
0 


pf i 2 sin(2k +1)aP(r)dz. 
Jo 


It should be noted that the set of functions 
V 2cos2ax, V2sin2axv, V 2sin4ar,--- 


which serves as a basis for the formation of the series (8.1), is merely the 
ordinary Fourier set transformed to the unit interval and normalized. 


II]. THE CASE p = 1 
Before proceeding to more general considerations we propose to discuss 
the series (8.1) and (8.2) when p 1 almost everywhere. In this case 
P(x) = x so that the expansions appear as series formed from the two 
normal orthogonal sets 


(9) 1, V2eos2ax2, V2sin2ax, V2cos4azx, 
(10) V2sinaz V2¢os3ax, V2sin3az,---, 


on the interval (0,1). Concerning the expansions of an arbitrary function 
integrable in the sense of Lebesgue we shall prove the 

THEOREM I. The formal expansions of an arbitrary summable function in 
terms of the sets (9) and (10) on the interval (0,1) are the formal expansions 
in terms of the normalized Fourier set for the interval (0,2) of two functions 
Ji and fe respectively. 

Proof. The normalized Fourier set for the interval (0,2) is seen to be 


1 
(11) 77s, cosax, sinaz, cos2nx, sin2axz, cos3az, sin4daz,- 


V2 


If we denote by ax, b; the Fourier coefficients of a summable function in 
terms of (11) we have 


41k 
By = 
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= = { f(x) + 


0 0 

| Jsin2k f(x) +f(2+1 )}sin2k 
e70 eJ0 


*1 
feos(2k+1)axrdx | f(a) —f(a +1)} cos(2k+1)aadz, 


) 


e7V0 


) 


If there is given an arbitrary summable function /(«°), 0 
define f,(.°) for the interval (0,2) by the identities 


f(x), 0: 
1) S(*), 0 


Then the Fourier series of /, on (0,2) has the coefficients 


1 1 a1 
lo V2 dex - 2 | box 2 | 
«70 «JV 


(ek+1 = 0. 


Furthermore the Fourier series thus set up, 


,< 
ax cos2kax- | Do sin2k ax}, 


is term for term the formal series for /() in terms of (9) on (0,1). 
Similarly, if we take /:() such that 
S(2), 


there results 


ao = Ax box = 


= 2 2k + 
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and the formal series for #2 becomes 


PAL -1cos(2k +1) art 
k=0 


which is identified at once as the formal expansion of (2) in terms of 
(10) on (0,1). 

CorOLLARY I. The term-by-term difference of the formal series for f(x) 
in terms of the Fourier set (9) and of the set (10) converges uniformly to 
zero on any closed interval interior to (0,1). 

Proof. This term-by-term difference is the Fourier series on (0,2) for 
since 0,0<-2z 1, it follows that this Fourier series 
converges uniformly to zero on any closed interval interior to (0,1), by one 
of the most elementary properties of Fourier series. 

CoroLLaryY II. Jf f(a) ts continuous and of bounded variation in a neigh- 
horhood on the right of x Q and is continuous and of bounded variation 
in a neighborhood on the left of x = 1, the formal series for f(x) in terms 
of (10) converges at x = 0 to 4(f(0+)—f(1—)) and at x = 1 to 
(—f(0+) + f—)). 

Proof. This is a direct consequence of the fact that the series is a Fourier 
series for fo(a) on (0,2). It is sufficient to recall that under the conditions 
stated this Fourier series converges at « = 0 to $(f(0+)+A4(2—)) 
and at 2 = 1 to 4(f2(1—)+A(1+4+)). 

It is clear that a large number of properties of Fourier series will be 
carried over by Theorem I and Corollary I to the series formed from (10).* 
In particular there exists no summable function other than a function 
identically zero almost everywhere on (0,1) such that its Fourier-like co- 
efficients with respect to (10) all vanish. Again we see by Corollary I that 
all properties of divergence and of both uniform and non-uniform conver- 
gence and summability are shared alike on any closed interval interior to 
(0,1) by the expansions formed from the Fourier set (9) and the set (10). 


*It is of interest to notice that analogous considerations for the Fourier series on 
(0,27) will lead to the series in terms of the sine set and of the cosine set on (0,7), 


9 
sin2ar,°°:; cosz, 


The term-by-term difference of the formal expansions of any summable function on (0, 7) 
in terms of these two sets converges uniformly to zero on any closed interval interior to 
(0,7). For functions summable with summable square this was proved by Walsh and 
Wiener, Journal of Mathematics and Physics of the Massachusetts Institute 
of Technology, vol. 1 (1922), pp. 103-122, especially pp. 115-116. 
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Gibbs’ phenomenon is also common to the two series in the neighborhood 
of any interior point of the unit interval for which it occurs for the series 
from (9). The behavior of an expansion in terms of (10) on the entire 
unit interval may be discussed by means of the Fourier series for /. on 
(0,2). The set (10) is thus of considerable interest as a set of trigonometric 
functions distinct from the Fourier set but resembling it very closely and 
arising from a similar differential system. We may also observe that ex- 
pansions in terms of (10) can be investigated separately from the Fourier 
series by means of an integral analogous to the Dirichlet integral which 
is central in the theory of the latter. 


III. THE TRANSFORMATION = P(x) 


It is evident at a glance that the transformation ¢ = P(2) is of fun- 
damental importance in the discussion of the series (8.1) and (8.2) since, 
no matter what summable function is used in calculating the coefficients, 
any sum attributed to the series must be of the form f[P(a)]. The 
relation of this sum function, if it exists, to the function for which the 
formal series is written down will depend on the nature of this transformation 
and its application to the integrals which are the coefficients of the series. 
We therefore devote a section to a detailed study of the relation ¢ = P(x). 

From the assumptions made in § I concerning the function p() we may 
immediately deduce the following theorem which describes the manner of 


x 
oscillation of the continuous function = paz: 


THEOREM II. The interval (0,1) can be divided into a finite set of non- 
overlapping intervals &m such that 


(a) the set covers the whole unit interval; 


(b) on & the function P(2) either always increases or always decreases and 
Surthermore P(x) lies between consecutive integers: k < P(x) << k+1, 
k = —M,—M-+1,---,—1,0,1,---, +H; 

(«) if on & the function P(x ) always increases (decreases) and lies between 
consecutive integers k and k+-1, then on the two adjacent &-intervals 
P(x) either always decreases (increases) or does not lie between the 
integers k and k-+1. 

Proof. On any one of the intervals 9; of § I, p(x) is either positive 
or negative almost everywhere. For the purpose of the argument we may 
assume the first possibility to be the case. Then for any two points 2;, 

*Lo 
on that interval we have = pdx >, 
since there exists on (7,, 7.) a measurable set of ‘points with measure 
greater than zero on which p(x) is positive. The same type of reasoning 
applies when p(x) is negative on g; rather than positive. 
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Each interval @;, on which P(x) is now known to be either always 
increasing or always decreasing, can be divided into a finite number of sub- 
intervals on each of which P(x) lies between consecutive integers and on 
no two of which P(a) lies between the same pair of integers. The number 
of sub-intervals thus determined is finite since P( 2) is continuous. 

The sub-intervals just defined form a set of intervals §,, &,---, &, with 
all the properties set forth in the statement of the theorem. 

It is now an exceedingly simple matter to determine the manner in which 
the intervals §; are transformed by the relation ¢ = P(x). We have 

THEOREM IIT. The transformation t = P(x) takes the set of intervals 
&; into a set of intervals % such that 

(a)t; has positive or negative sense according as P(x) increases or de- 

creases on §&;; 

(b) the algebraic sum of the intervals 1; is the unit t-interval; 

(c) tj lies between successive integers. 

Proof. The interval £;:(a;,;) clearly goes over into the interval 7;: 
[P(a), P(bi)]. Statements (a) and (c) are then made apparent by 
Theorem II (b). It is also easy to see that if 0, 7, v2,---,%,~—1,1 are 
the end points of the intervals g; of § I the algebraic sum of the intervals 7; is 
found by adding the intervals [P(0), P(m)], [P(m), P(z2)], 
[P(an-1), P(1)] algebraically. The sum is then seen to be the interval 
[P(0), P(1)] which is precisely the interval (0,1). 

Before going on to discuss the transformation of integrals, we shall state 
two necessary lemmas of a general nature. Of these we shall prove the 


Lemma I, If ¢ = t(x) és the Lebesgue integral a’ of a 
summable function g(x) which is positive almost everywhere a < x < b, then 
(a)t = t(ax) has a unique inverse x = x(t) which is a one-valued always 
increasing function of <t< 
(b)7(t) is equal to the Lebesgue integral 


a’ 


If t = t(x) is the Lebesgue integral . g(x)da+a’ of a summable func- 


tion g(x) which is negative almost everywhere a < x < b, then 


(a)t = t(x) has a unique inverse x = x(t) which is a one-valued always 
decreasing function of t, b) <t< a’; 
(b) x(t) is equal to the Lebesgue integral J 


at 


a 

Proof. As we have already seen in the proof of Theorem II, ¢(z) is 
an always increasing continuous function of x when g(x) is positive almost 
everywhere on (a, )); statement (a) of the first part of the lemma follows 
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immediately from this fact. Turning now to (0) in the first part of the 
lemma, we examine the equality 


t+h )—zx(t) = h’ 


h t(a+h')—t(x)’ 
where h,h’ vanish together. The reciprocal of the term on the right 
has the limit g(x) almost everywhere; since g(x) is positive almost every- 
where the term itself has a finite limit almost everywhere. Corresponding 
to the set of zero measure on the x-axis for which the limit does not exist 
or is not equal to g(x) there is a set of zero measure on the ¢-axis. To 
establish this statement we note that the measure of the set of points on 
the ¢-axis corresponding to a finite or denumerably infinite set of non- 


overlapping intervals E on the z-axis is equal to the integral J pac. 
E 


We can now enclose any set of zero measure on the z-axis in an open 
point set of measure less than 6*, that is, in a finite or denumerably 
infinite set of non-overlapping closed intervals E+. By choosing 6 sufficiently 
small we can make the measure of £’, which corresponds to E and 
therefore contains the set corresponding to the set of zero measure on the 
x-axis, less than any preassigned positive «: 


m( = [ m(E)< 0. 

JE 
Thus we must have the result that lim,_,,.(#(¢+h)—2(t))/h exists and 
is equal to 1/g[x(t)] almost everywhere, a’ < ¢ < b’. 

If we can show that z(t) is a Lebesgue integral we can then state 
that the function is the indefinite integral of its differential coefficient, 
c= | — +a. It is well known that a necessary and sufficient 

Ja’ glx(t)] 
condition that a function be an integral in the sense of Lebesgue is that 
it be absolutely continuous. We shall prove that this is true for x(t); 
that is, we shall show that >| x«(b;)—2(a;)| formed for any set of non- 
overlapping closed intervals E’ can be made uniformly small by taking 
m(E’) sufficiently small. In other words, if we are given any positive «, 
however small, we can find a positive 6 such that if m(#’)<o then 
m(#)<«, where £ is the set of intervals on the z-axis corresponding 
to EZ’. For the set of points on which g < Vd is of measure 1(¢) 


vanishing with 6; we pick VY and require m(E’) <4. 


*de la Vallée Poussin, loc. cit., § 21. 

7 de la Vallée Poussin, loc. cit., § 14. 

t Vitali, Atti della Reale Accademia di Torino, vol. 40 (1905), pp. 1021-1034; 
Lebesgue, loc. cit., p. 129, footnote. 


1924] AN EXPANSION PROBLEM 


Splitting E into the measurable sets of points £,, Z, on which g> 
g V 6 respectively, we see that 


Vd m(E;) < | gdx | gdx = m(E")<4, 
JE, JE 


whence 
m( EF) = m(F;)+ m( £2) 


The proof of the first part of the lemma is thus completed. The second 
part, in which g(a) is negative almost everywhere, is proved in the same 
manner. 

The second lemma we state without proof. It is 

LemMA Ifa = x(t) is the Lebesgue integral of a summable function h(t) 
which is positive (negative) almost everywhere on (a’, b'), then for any sum- 
mable Junction I(x) 


fa) de ‘(t)] - da = h(t) dt, 
b = 2z(V’), a = x(a’). 


With the aid of Lemmas I and II we are now able to demonstrate 
THEOREM IV, Jf p(x) f(x) is summable, 0O< a4 <1, then 


where X;(t) is the inverse of t= P(x) for the interval &;; and, conversely, 
if f(t) is summable and is defined outside the unit interval by the periodic 
relation f(t +1) = f(t) or by the anti-periodic relation f(t+1)= —/(t), 
then 


da = ) de. 


JT; e 

Proof. On the interval §;, the function ¢ — P(z) is the Lebesgue integral 
of a funetion positive (negative) almost everywhere on that interval, so 
that by Lemma I we can invert t= P(x) by means of the function 


at 

dt 
+ 


* Lebesgue, Annales de la Faculté des Sciences de Toulouse, ser. 3, vol. 1 (1909), 
pp. 25-117, especially p. 44; Hobson, The Theory of Functions of a Real Variable, second 
edition, Cambridge University Press, 1921, § 440, pp. 592-595. 
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for ¢ on 7%: (aj, bj) and x on §&:(a, b;). We then apply the second Lemma 
to each of the two integrals | praz, |_ fat remembering that 1/p[X;(t)] 
is positive almost everywhere or negative almost everywhere on 7; according 
as p is positive or negative almost everywhere on §;. It should be noted 
that since p is bounded, p/ is certainly summable if / is. 

In order to pass to the consideration of integrals like those which appear 
as the coefficients in (8.1) and (8.2), we introduce new concepts. If we 
define X;(t) by the periodic relation X;(¢-++ 1) = Xi(t) for all the intervals 
ti+k, k=O, +1,--- then for some k, say k;, one of these intervals will 
lie entirely on the unit interval and thus 


1X00) dt = | at. 


The new interval 7;-+ 4; will still have definite positive or negative sense; 
to eliminate this circumstance, we let 7; be the interval 7;-+ 4; taken in the 
positive sense. We can then write 


dt = [sen f[Xi(t)] dt 


e i 
where sgn 7; = +1 if t% has positive sense and sgn 7; = —1 if 7; has 
negative sense. With a view to the discussion of the series (8.1) we now 
define a function 


gi (t) Jiu (t) + fi2(t) + + fim (t), 


where fi; = sgn 1; [Xi] on tj and is identically zero elsewhere on (0, 1). 
Similarly, with a view to the treatment of (8.2) we define a second function 


= f(t) + foo(t) + +fom(t), 


where fo; = (—)* sen x; f[Xi] on ct; and is identically zero elsewhere 
on (0,1). These two functions are so related to /(«) that we shall call 
them the first and second P-average functions for f(x), respectively. The 
introduction of these functions is prompted by considerations which will be 
more apparent later in our work. 

We now prove two theorems on integrals which will serve as the basis 
for our analysis of the series (8.1) and (8.2). We begin with 
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THEOREM V. If p(x) f(x) is summable, then g,(t) and gs(t) are sum- 
mable and 


*1 
f(2) dx = | g(t) dt. 
0 e 


If p(x) f?(x) is also summable, then g,(t) and gs(t) are both summable 
with summable square. If f(x) is bounded, then g,(t) and gz(t) are also 
bounded. 

Proof. If p(x) f(x) is summable, then, by Theorem IV, 


*1 
= | Su(t) dt 


*1 


| pfdx = (—1)k | Su dt 
J 


so that fi; and f; are summable. Since gi (¢) is a linear combination of 
the fic and g2(t) of the both gi(t) and g(t) are summable. The 
addition of the equalities 


ps dz =f Su dt 
E; 0 
41 1 
J = dt. 
0 0 


If pf? is summable also, then Jor dx exists and by Theorem IV and 


gives us the fact that 


the definition of the functions /;, fo: is equal to 


41 
| f2, at dt. 
eJ0 0 


Thus f2 and /? are summable so that this is also true of g, and gs. 


1 
We may note that if f pf? dx exists the inequality 


1 | 41 
| psf dz <{ dx 
0 
Vip| Vip! dx <| | | lp| dz 
0 J0 0 


1 
proves the existence of f pfdz. 


25 
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Lastly, if # is bounded so are /;; and f2;. In consequence g, and g, 
are also bounded. 

The other theorem serves as a sort of converse io the one we have 
just proved. Its statement follows: 

THEOREM VI. Jf h(t) is summable and is defined outside the unit interval 
by the periodic relation h(t+1) = h(t) or the anti-periodic relation 
h(t+1) = —h(t), then p(x) h[ P(x)] és summable and 


+1 
ph(P| dz = | h(t) dt. 


0 0 


If h(t) is in addition of summable square, then ph*| P| is summable. If h(t) 
is bounded, then so is h[ P(x)]. 
Proof. In any of the three cases we can write 


m m 


h(t)dt= h(t)dt >| p(a)h[P(x)] dx dz. 
1 Ti E; 


le 


Thus ph| P] is summable. In case we have /?(¢t) summable we may replace /(t) 


by in the equalities just written down, obtaining /°(t) dt = | 
of 


and in this way proving that ph*|P] is summable. It is obvious that if (1) 
is bounded, h[{P] is also bounded. 

With Theorems V and VI we are prepared to treat the two series (8.1) 
and (8.2). It seems in place, however, to say a word about the possibility 
of extending the results of this section to transformations associated with 
a somewhat more general class of functions p. To obtain analogues of 
Theorems II and III it is sufficient to require that it be possible to deter- 
mine a denumerably infinite set of non-overlapping intervals 0, 02, @s, - 
covering the unit interval except for a set X of zero measure and such that 
on @; the function p() is either positive or negative almost everywhere. 
The analysis necessary to establish a theorem parallel to Theorem III is 
somewhat more complicated than that which we have found sufficient for 
the type of function p(s) we are considering. Theorem 1V remains the 
same. By investigating some simple convergence questions we can set up 
the P-average functions for any function /(), but Theorems V and VI are 
no longer true, at least for the general class of functions p(a’) just defined. 
To bring out the essential points of the treatment which we have outlined, 
we shall take up a simple example. 

We shall define p(x) so as to be constant on each of a denumerably 
infinite set of intervals @,, @:,---. It is well known that the infinite series 
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1+2>'%(1/n*) converges to a positive constant C. The interval @, shall 
be of length 1/C’ and shall have « 1 as an end point. Then g, shall 
abut on @, and be of length 1/4’, gy shall abut on g, and be of length 1/4. 
In general, @::1 shall abut on @;, and @%, @x41 shall each be of length 
i((k+1)?C). Then S* m(e) 1. We shall denote the end points of 


the intervals @,, @.,--- by 1 ly >> >> NOW we put 


p(x) +C, 
p(2) 


p (ax) 0 


A few simple computations show that P() - { p(«)d~ is a linear function 
) 


of on each g-interval, while P(ax) = Wk +1)2, k=0,1,2,--- and 
P(dxii) = 0, k = 0,1,2,---. The o-intervals of this example are 
precisely the §-intervals analogous to those of Theorem II. By the trans- 
formation ¢ = P(x) there corresponds to the interval g24; the interval 
t+1:(0, 1/(k-+1)*) taken in the positive sense, and to the interval ox the 
interval t2,:(0, 1/(-+1)*) taken in the negative sense. Since all the inter- 
vals 7; are on the unit interval, the two P-average functions coincide in 
this case. In constructing this one we will obtain instead of a finite sum 
an infinite series which, however, has the peculiarity that on any interval 
(e, 1), ¢> 0, it is actually a finite sum. As an example we may take the 
bounded summable function which is identical to 1 on @,, 03, @5,--+, and 
identical to zero elsewhere. We have at once 


| pdz = | dt. 


Thus if we put fox41(t) 1 on f,41 and let it vanish identically elsewhere 
we find that the P-average functions for this particular /() are given by 
the infinite series g(t) oS (4). It is quickly seen that g(t) =k 


on [1/(k Now f(x), being bounded, is of summable square; 
but, on the other hand, we see that 


41 D 
2 > 2 


a divergent series. Thus g(t) is not of summable square, and Theorem V 
is not true as a whole for this transformation. In the same way we can 


25* 
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1 1 
show that, although 4¢ * is summable, }p(a)P * (x) is not; for if 
1 


i» P * were summable its integral over the set of intervals @1, 03, 05,--. 
1 


oo 
would exist and would be given by » f pP * dx. This last expression 


ON 


turns out to be a divergent series since 


Thus Theorem VI is not true as a whole for this special case. 

The example which we have just discussed makes it clear why we have 
not attempted to take p(s) of more general character, and serves to 
indicate along what lines any generalization can be expected to take place. 
With this we may leave the developments of this section. 


IV. THE BEHAVIOR OF THE FORMAL EXPANSIONS 


Instead of discussing the series (8.1) and (8.2) alone we shall consider 
a wider class of series including these two as special cases. If we are 
given any normal orthogonal set gi,(¢), where gi.(t+1) = g(t), we 
can apply Theorem VI to the summable function y(t) gix(t), whence 


‘1 itk 
j guigudt = P Pil dx ={ 
e/J0 J 0 1 i=k 


The set (9) is a special case. Similarly, if we are given a set of normal 
orthogonal functions where go.(t-+1) = — yx (t) we find 


Jo? Pl yal Plde 


The set (10) is a special example. It is natural to form from the sets 


gulP], gis|P], 
yal P), y22[P], 


formal series 


41 
of the class for which ay, = DS PI dx and ax, = i} pS Gx | dz 


az P| and ax respectively for any function 
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exist. If gix, 92, are bounded then this class of functions is at least the 

class of all functions for which p(a)f(2) is summable. If gx, Gx are 

summable with summable square then this class is at least the class of all 

functions for which p (x) f(z) is summable. By Theorem VI the Lebesgue 
1 


integrals p¥,[P] dx exist. Because of the existence of 


“1 
| the inequalities 


| 
DA < PS | P| 


1 


| 41 
PI P| dx) < DA 


a1 
p go. | ics p 95, da 


demonstrate the summability of the measurable functions pfgu, pf px. 

We come next to 

THEOREM VII. The Fourier-like coefficients of f(x) in terms of the set 
of functions gu[P], ¢12[P),--- are the Fourier-like coefficients of the first 
P-average function for f(x) in terms of the set 9u(t), gi2(t),--+3 and the 
Fourier-like coefficients of f(a) in terms of the set gu[P], 92[P],--- are 
the Fourier-like coefficients of the second P-average function for f(x) in 
terms of the set goi(t), 

Proof. If Gu, (#t) is the first P-average function for the product f9%[P] 
then, by Theorem V, Gy.dt psoul P\dt. Now when evalu- 
ated is precisely gi (t)¢ix(t) where g; is the first P-average function for 

1 
f(x). Hence ay, = pf oul P\dz = (t)ou(t)dt. In the same way, 
if G(t) is the first P-average function for the product fyx[P] and g2 is 
the second P-average function for f(a) then 


1 


Go(t) = go(t) pox(t), = PS = | 92 (t) pox (t)dt. 
Jo ‘ 


CoroLLARY I. Jf pf? dx exists, a series formed from P| by 
e 
grouping its terms suitably will converge essentially uniformly to gi| P| where 


y 

i] 

0 
| V 
J0 

e 

al 

ts 

on 

ly 
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gi(t+1) = gi(t); and a series formed from dx gor | P| by grouping 
its terms properly will converge essentially uniformly to g2[(P] where 
gz (t+1) = 

Proof. By the hypothesis that f pf*dx exists we have the fact that 
the two P-average functions for /() are summable with summable square, 
Consider the first of these, g(t). By the Riesz-Fischer theorem, in the 
form given to it by Weyl, there exists a series formed by grouping the 
terms of the series >° aix91%(t) which will converge essentially uniformly 
to a function gi(t), where g:(¢)—gi(t) is a function all of whose coef- 
ficients with respect to yu(t), gi2(t),--- vanish.t By essentially uniform 
convergence we mean that, given any decreasing sequence of positive ¢’s 
with limit zero, ¢, > 6 > 6 >---, we can find on the unit interval a sequence 
of measurable point sets Ey, Hy, Ey,--- each contained in the succeeding 
set, with m( Ey) > 1—e;, such that on each set of the sequence the series 
converges uniformly to its limit. On the unit z-interval we now have a 
sequence of point sets L,, Ey, E;,--- where E; is the set of points for 
which ¢ = P(x) takes on values congruent (modulo 1) to a point of Ey. 
We see that on E; the series formed by grouping the terms of SY? ax ¢1%[ P] 
converges uniformly to gi[P] because of the periodicity of the functions 
gu(t). Now CE, the set complementary to E;, consists of all the points 
of the unit «-interval for which ¢ — P() is not congruent (modulo 1) to 
a point of Ey; that is, of all the points for which ¢ = P(x) is congruent 
to a point of CE/. Thus CE, can be obtained as the sum of the point 
sets @2k, @3k,°**, Where is the set on &; corresponding by the 
relation P(a) to a set congruent (modulo 1) to CE. Sincem(CEY) < 
can enclose CE; in an open set Oj of measure < 2¢,. On the interval &; the 
set oj, corresponding to O; contains ex. Since for values of ¢ on 7% we have 


_dt 


it follows that m(ox) is given by the Lebesgue integral of sgn 7;/p[ Xi] 
taken over the points of O; on 7}. From this point on we apply the 
arguments we used in discussing a similar situation in the proof of Lemma I. 
By choosing an ¢ small enough we can make m(CE,) = >” ,m(ei) 

- >'_, m(oix) less than any preassigned positive «. Thus, given a sequence 


of decreasing positive ¢,’s with limit zero, we can find a subsequence of 
E,, E2,---, say Ey’, Ey’, Ej’,---, such that m( Ey’) >1—e. The 


7 Plancherel, Rendiconti del Circolo Matematico di Palermo, vol. 30 (1910), 
pp. 289-335, chap. I. 
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series formed by grouping the terms of >*° auxg1%[{P] therefore converges 
essentially uniformly to the function gf[P] where gi(¢+1) = gi(t). If 
the set y%(¢) is closed with respect to all functions summable with sum- 
mable square—that is, if the only function of this class which has all zero 
coefficients with respect to the set is identically zero except on a set of 
measure zero—then gi(t) = gi(t) almost everywhere. Exactly the same 
sort of reasoning applies to the function and the series ax g(t), 
the only difference being that here the functions y(t) are anti-periodic 
instead of periodic. 

CoroLuary Il. Jf the set gix(t) is closed with respect to all functions 
summable with summable square, and if fi(a) and fo(x) are two functions 
such that “Ph dx and pAzde exist, then a necessary and sufficient 
condition that fila) and fo(a) have the same coefficients with respect to 
gulP], P|,--- ts that their first P-average functions coincide almost 
everywhere. If the set pix is closed with respect to all summable functions, 


and if fila) and fo(x) are two functions such that phi dx and phodsx 
exist, then a necessary and sufficient condition that fi(x) and fo(x) have 
the same coefficients with respect to the set gu[P], ¢i2[P],--- is that their 
first P-average functions coincide almost everywhere. Similar statements 
involving the second P-average functions can be made for the set 9u[P), 
PI, 

As we have already remarked in § II the sets (9) and (10) are closed 
with respect to all summable functions. 

Corollary I to Theorem VII suggests a theorem which we shall prove 
independently as 

THEOREM VIII. Jf the function P(x) is not monotone, there are infinitely 
many non-null functions whose P-average functions are both identically zero. 

Proof. By a non-null function we mean any function which is not a null 
function; and a null function is any function identically zero except on 
a set of measure zero. When we assume that P(x) is not monotone we 
assume that there are two or more of the intervals og; described in § I. 
For the sake of definiteness we shall assume that on @; : (0, a) the function 
P(x) inereases from zero to P(a) while on @2:(a, 6) it decreases from 
P(a) to P(b). If P(b) = 0 we consider the interval A; made up of @2 and 
the sub-interval of 9: for which P(b) < P(x) < P(a); in case P(b) < 0 
we take the interval A. composed of g; and the sub-interval of @2 for which 
0< P(x) < P(a). If there is given any bounded summable function 
h(t) defined for all values of ¢ on the interval 0 < ¢ < P(a), we set 
J(«) identically equal to h[P] on Ri or Ry as the case may be, and 


identically equal to zero elsewhere, 0 < « < 1. It is then not difficult 
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to see by reference to the definitions of the two P-average functions for 
f(x) that both these functions vanish in the case of the bounded summable 
function which we have constructed. All the coefficients of f(x) with 
respect to the sets gu[P], gi2[P],--- and ga[P], 22[P],--- exist and 
vanish by Theorem VII. It is, of course, unnecessary but simpler to confine 
ourselves to bounded functions. 

Theorem VII, Corollary I, and Theorem VIII reveal the striking peculiarities 
of the expansions of §I and their generalizations; peculiarities entirely 
new, so far as we are informed, in expansion problems. For any function 
f(x) such that p(x) f?(x) dx exists, the formal series will always represent 
in a certain sense a definite function, which will be of the specialized 
form gi[P] or g2[P]. It is clear that if for dx exists then f(a) can 
be represented by its formal series in terms of a closed set gu:[P], 912[ P], °°: 
if and only if its functional values are so distributed that it is of the form 
mlP], m(t+1) = m(t) except possibly on a set of measure zero; and 
by its formal series in terms of a closed set ga[P], g22[P],--- if and 
only if its functional values are so distributed that it is of the form 
golP], go(t+1) = —ge(t) except possibly on a set of measure zero. 
There are infinitely many functions all having the same formal expansions 
in terms of the two general sets of functions which we are discussing, 
and this quite apart from any considerations of the closure of the sets 
u(t), pio(t),--- and go(t), go(t),---. The remarks which have just 
been made are, of course, expected to apply only in case P(x) is not 
monotone. If P(s) is monotone the expansions have no unusual properties. 

The detailed consideration of a formal series here is thrown back on 
to that of the series for gi(t) or the series (t) 
for g2(t). We do not intend to go deeply into any such discussion, but 
will close with one theorem bearing on the series (8.1) and (8.2): 

THEOREM IX. If the function f(x) is of bounded variation and if its 
first P-average function is extended by the periodic relation g:(t+1) = g(t) 
then the formal series (8.1) converges at every point x for which P(x) = to 
to the value i(g(t.+0)+m:(t—O0)); and if its second P-average 
Junction is extended by the relation g.(t+1) = —ge2(t) then the formal 
series (8.2) converges at every point x for which P(x) = ty to the value 
2(g2(to +0) + g2(tr-—0)). 

Proof. The function f(z) can be represented as the difference of two 
monotone increasing functions f, and f2 because it is of bounded variation. 
Now the functions /i;(¢) and fo; (¢) used in defining the P-average functions 
for f(x) are constant multiples of /[X;(t)] on cj and zero elsewhere. X;(t) 
is periodic of period 1 and is monotone on z;. Since a monotone function 
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of a monotone function is monotone, both f,[Y;(t)] and #2LX; (t)] are 
monotone on and f[X;] = is therefore of bounded 
variation on 7]. Thus g,(t) and gs(t) are also of bounded variation, being 
jinear combinations of functions of bounded variation. Now the series for g, 
in terms of the set (9) is the Fourier series for a function of bounded variation 
and therefore converges to $(g; +0)+ 9:(t—0)) at The substitution 
of P(x) for ¢ gives the first part of the theorem as stated. Likewise, as 
we saw in Theorem I, the series for g.(t), where g.(t+1) = — ge(t), in 
terms of (10) is the Fourier series on the interval (0,2) for a function 
of bounded variation. The series converges at the point ¢ to the value 
+0) + ge(t—O0)). The substitution of P(x) for ¢ completes the 
theorem. 

Other theorems of a similar nature will readily suggest themselves to the 
reader. We shall not go into any further detail, having obtained by the 
work up to this point a fairly complete idea of the character of the expansion 
problems arising from the two differential systems (1) and (2). 

HARVARD UNIVERSITY, 

CAMBRIDGE, MAss. 


ON THE INDEPENDENCE OF 
PRINCIPAL MINORS OF DETERMINANTS* 


BY 


E. B. STOUFFER 


INTRODUCTION 

A principal minor, or coaxial minor, of a determinant A is a minor 
obtained by striking out from A the same rows as columns. There are 
2”—1 principal minors of a determinant of the nth order, the determinant 
itself being included, but only n?—n-+1 of them are independent.+ 

For n = 1, 2, 3 the principal minors are all independent and for n = 4 
the relations between them have been quite extensively studied.t For 
n> 4 little has been published either as to which minors constitute an 
independent set or as to the relations between the minors. 

It is one of the purposes of this paper to determine several different 
types of complete sets of independent principal minors of the general 
determinant of the nth order and to show how the elements of the determinant 
may be expressed in terms of the minors of an independent set. 

If we have a second determinant of the mth order with elements 
independent of the elements of the first determinant, there is a definite 
set of determinants obtained by replacing one or more columns of the first 
determinant by the corresponding column or columns of the second deter- 
minant. The set of principal minors of all such determinants is greatly 
enlarged over the set from the single determinant. It is the second purpose 
of this paper to determine several types of complete sets of independent 
principal minors of this enlarged set. There is also determined in this 
paper a complete set of independent principal minors of the determinants 
obtained when the above process is extended by adjoining to the original 
determinant more than one determinant of the nth order with independent 
elements. 


* Presented to the Society, April 14, 1922, and December 1, 1923. 

+See MacMahon, Philosophical Transactions of the Royal Society of 
London, vol. 185, pp. 111-160; Muir, Philosophical Magazine and Journal 
of Science, ser. 5, vol. 38, pp. 537-541. 

t See, e. g., MacMahon, loc. cit.; Nansen, Philosophical Magazine and Journal 
of Science, ser. 5, vol. 44, pp. 362-67; Muir, Transactions of the Royal Society 
of Edinburgh, vol. 39, pp. 323-339. 
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It is the third purpose of this paper to prove the independence of certain 
sums of principal minors. In the fourth section of the paper the independence 
of these sums is used to determine the possibility of expressing polynomials 
as determinants with linear elements. 


I 
Let A represent a determinant of the mth order in which the element 
in row 2 and column 7 is denoted by a. Furthermore, let (Agim...) denote 
the principal minor obtained by striking out in A all columns and rows 
except those numbered k,/, m,n,---. Thus,* 
| Ay | 


(Ai) = au, (4y) = | 


| Aji Ajj | 


It has been shown by Muir? that the values of the principal minors of 
a determinant A with general elements aj are not affected if certain sets 
of »—1 of the elements are put equal to unity. For instance, if we 
divide the elements of the second column by as, of the third column by a,;, 
of the fourth column by a,,, and so forth, and also multiply the elements 
of the second row by as, of the third row by as, of the fourth row 
by a4, and so forth, there is no effect either upon the values of the 
principal minors or upon the independence of the individual elements, 


except that aj = 1(y = 2, 3,---,m). The fact that there are only 
n?>—n-+1 independent elements remaining is Muir’s proof that there cannot 
be more than n?—n-+1 independent principal minors. 

It may be shown similarly that there is no loss of generality in assuming 
ayy = 1 in place of ay = 1. 

THEOREM I. The minors 
(1) (Aj), (Ay), Gj = 1,2,---,n;2 <j) 


constitute a total of n®—n-+1 independent principal minors. 

We shall prove this theorem by showing that the functional matrix of 
these minors with respect to the elements aj contains a non-vanishing 
determinant of order n*—n-+1. In order to simplify the process we 
shall assume that ay; = 1 (j = 2,3,---,) and also that, after the partial 
derivatives with respect to the remaining aj have been formed, the elements 
aij > 7 > 1) are all put equal to zero. Then all an: except 


* We shall assume throughout this paper that the subscripts in the symbols for principal 
minors have been arranged in order of ascending magnitude. 
7 Muir, Philosophical Magazine and Journal of Science, ser.5, vol. 38, 


pp. 537-541. 
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that (A;)/8a% = 1; all 6(Ay)/@am = 0 (1<m,m4l) except that 
0(Ay)/8 ag = —ay; all 6(A1y)/Oam = O(1 <m <1) except that 
0 (A4)/9 ai; = aj1. Under these conditions there is in the functional 
matrix a determinant of order n?—n-+-1 which may be expressed as the 
product of the three determinants 


(2) — 


[(Ai2), ( Ais): (Ain), ( Ags), (Aan), ( An—1,n)] 


Aen, 35, ***s 


These determinants are all different from zero since in each of them the 
elements of the principal diagonal are different from zero and all other 
elements are equal to zero. The minors (1) are therefore independent. 

The explicit expressions for A and for the principal minors of A in terms 
of the minors (1) are obtained by expressing the elements of A in terms 
of the minors (1). Under the condition that a; 41 (j >1), we have 


(4;) = au, 
(Ayj) = a1 aj, 
(Ay) = (4>1), 
= au Ay— da aj + aj Gj + da aj 


whence 
au = (A,), 


ay a = (A;) (Aj) — (Ay) (2 >1), 
Aj + din = (Avy) —(A1) (Ay) + 2 (Ad) (Aa) — (dj) (Ara) — (Ad) 


(6) 


The expressions for all the elements of A in terms of (A;), (Ay), (Ary) 
follow immediately from these equations.* 

Another complete set of independent principal minors of A is obtained 
if the minors (Ay;) of the set (1) are replaced by the minors (A19;), (A12a;), 


* MacMahon, (Philosophical Transactions of the Royal Society of London, 
vol. 185, p. 147) states that it is impossible to express a determinant of odd order in terms 
of its principal minors. By the above process this statement is shown to be incorrect. 


(4) ), (Aisa), -, (Ayn), (Aisa), (Ais), (Aran), » (A1,n—1,n)] 
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(Ayssj), «++, (Args...n), Where 7 in each case takes all the values up to 
which are greater than the subscript immediately preceding it. The functional 
matrix will again be used in proving the independence of this set of n*—n-+1 
principal minors. If we assume as before that aj = 1 (j = 2,3,---, m) 
and that aj (i >j > 1) are put equal to zero after the partial derivatives 
are formed, the matrix contains a determinant of order n*—n-+1 which 
consists of the product of the determinants (2) and (3) and the determinant 


(7) 0 [( A128), (Aizs), » (Aizen ), ( Ajzas), | Aj23...n)] 


‘ a 
O[ aes, dos, +++, Aon, Anan] 


Under the above conditions 


1 1 1 1 
0 ( 
= tay Ags a3,i-1 asi + 0 (1 
0 Ajj | 
0 0 ai-1,i-1 Q-1,i 
and 
0 (m<1l;m>iorl>j). 
0 Ami 


Consequently, in the determinant (7) all elements above the principal diagonal 
are equal to zero and all elements of the principal diagonal are different 
from zero. We thus have 

THEOREM I], The minors 


(5) (Aj), (Ay), ( Ajyj), (Aj93;), 


form a complete set of independent principal minors of the general deter- 
minant A. 

It is easily verified that the process of obtaining the expressions for the 
elements of A in terms of this second set of principal minors involves 
merely the successive solution of sets of two equations, one of which is 
linear and the other quadratic. However, the pairs of equations to be 
solved are not so easily written as before, since each pair has a different 
form from those which precede it. 

Still another set of independent principal minors is obtained if the 
}(n—1)(m—2) minors (Aj) of the set (1) are replaced by the minors 
(Ares), (Agsa), (Asis), (Aosas), (Args---n)- In order to 
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prove them independent, let us assume that aj = 1 (¢-+1 =~), instead 
of aj = 1, and that aj(¢4 >j-+1) are put equal to zero after the partial 
derivatives in the functional matrix are formed. Under these conditions 
the matrix contains a determinant of order n?—n-+1 which consists of 
the product of two non-vanishing determinants of the same form as (2) 
and (3) and the determinant 


al (124, An—2, M4; 125 Ain ] 
The latter determinant obviously has all elements of the principal 
diagonal different from zero and all elements above the principal diagonal 
equal to zero. We thus have 

THEOREM II]. The minors 
(10) (As), (Ay), = 2, 3, 7 <n) 
form a complete set of independent principal minors of the general deter- 
minant A. 

In order to express the elements of A in terms of the minors (10) it is 
again necessary to solve rather complicated pairs of equations, each pair 
consisting of a linear and a quadratic equation. 


Let B denote a second determinant of the nth order with elements by. 
A series of determinants of the mth order may be obtained by replacing 
one or more columns of A by the corresponding column or columns of B. 
Each of these new determinants gives a set of 2"—1 principal minors. 

Let us indicate by (Ajj... Byst...) the principal minor which belongs to the 
above set and which contains elements from columns j,k, --- (i<j <k--+) 
of the determinant A and from columns 7, s, t,--- (r <s <¢t---) of the 
determinant B. It is evident that the numbers 7, 7, k,---,7r,s, ¢,--+ must 
all be different. 

Just as in the previous case, we can make n—1 of the elements, say 
ayj(j = 2, 8,-++,m) or equal to unity without affecting 
the values of the principal minors. It follows that there are not more 
than 2n*—n-+1 independent principal minors. 

THEOREM IV. The minors 


(11) (Aj), (Bi), (Ay), (Ay Bi), (Bu), (Aj), (Ay B;), (Ay Bi), ( Buy) 
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form a complete set of independent principal minors of the determinants 
obtained by combining A and B. 

In order to prove the set (11) independent let us first observe that an 
element @im, bim (1 m) can occur only in minors having a subscript equal 
to 7 and another subscript equal to m. It follows that the functional 
matrix of the set (11) contains a determinant of order 22?—n-+1 con- 
sisting of the product of the determinants 


) (A,)] 


(12) 
Ann | 


and of determinants of the forms 


[(Au), (ArBi),(Bu)] Arg), (Ar By), (Bul 
bi | ( Sts yj). 


(13) (i, bu, bit a dij, ji, bi, 
The determinants (12) are already known to be different from zero. If, 
after forming the determinants (13), we put aj—0O(i>j>1) and 
b; =O(4 > 7), it is immediately evident that each of them is equal to 
the product of the elements in its principal diagonal, all of which are 
different from zero. 

By analogous processes we may prove 

THEOREM V. The set of principal minors 


(4;), (Bj), (Ay), (Ay Bj), (Bj), 
(14) 


1,i,j) Ay: (Aj, 2,3,..-, i-1 By), 
(i > 2), 


and the set 
(15) (Aj), (Bj), (Ajj), (4; Ber), 


(r = 2,3,---,n—1;itr<n) 


each form a complete set of independent principal minors of the determinants 
obtained by combining A and B. 

It should be noted that the sets (11), (14) and (15) are obtained 
from the sets (1), (8) and (10), respectively, by adjoining to each set 
the minors obtained by replacing in the determinants of the set the elements 
of the last column, the elements of the last two columns, and the elements 
of the first column and the last two columns, by the corresponding elements 


0[(A,). (A. O[(B,), (Be), (By)] 
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bj, except that only a part of the minors of the second order thus obtained 
are used. 
The n? minors 


(16) (Bi), (A; ( By) = 23 < 


when arranged in the above order each contain an element bj not in the 
minors which precede it and not in the minors of A. We thus have 

THEOREM VI. The minors (16) together with any complete and independent 
set of principal minors of A form a complete set of independent principal 
minors of the determinants formed by combining A and B. 

The process of expressing the elements aj and bj in terms of the minors 
of any one of the above independent sets involves the simultaneous solution 
of linear and quadratic equations. The process can be carried out quite 
easily for the set composed of the minors (16) and (4A;), (Ag), (Aww). 

The above process of forming determinants by combining two deter- 
minants A and B may be extended by combining /& determinants A, B, 
C, D,--- of the mth order. Select for the first column of the new 
determinant the first column of one of the k determinants, for the second 
column the second column of one of the k determinants, and so forth. 
The principal minors of all such determinants of the nth order will now be 
considered for the purpose of determining a complete and independent set. 

The symbol (Aix... Byst... Cuvw-...) Will be used to indicate the principal 
minor containing elements from columns 7, 7, k,--- of A, from columns 
ry, 8, t,--- of B, from columns u, v, w,--- of C, and so forth. 

Since »—1 of the kn? elements may again be assumed equal to unity, 
there cannot be more than kn?—n-+1 independent principal minors 
obtained from the whole set of combined determinants. 

Among the minors of the set (11) there is a group in which elements 
of B appear in the last column only and another group in which elements 
of B appear in the last two columns only. If in each of the minors of 
the first group the elements of the last column are replaced by the 
corresponding elements of C', D,---, successively, and if in each of the 
minors of the second group the elements of the next to the last column 
are replaced by the corresponding elements of C’, D,---, successively, we 
have, for each determinant C, D,---, n? new principal minors. That 
each minor thus found is independent of the others and of the minors of 
the set (11) is evident from the fact that, with the proper arrangement, 
each of them contains an element cj, dj,--- not in any of the preceding 
minors. Moreover, since the number of independent principal minors thus 
obtained is kn?—n-+1, the set is complete. 
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By exactly the same process we can obtain complete sets of minors for 
the enlarged system from either (14) or (16). The set (15) yields a 
complete set for A, B, C, D,--- if the elements of the last column of 
the minors which contain a single column of elements from £& and the 
elements of the first column of the minors which contain two columns of 
elements from B are replaced by the corresponding elements of C, D,---, 
successively. 

The minors to be added to any one of the sets (11), (14), (15), (16) 
may take many other forms, but the above systems are especially simple 
because the elements of C, D,--- appear only linearly. The fact that 
all elements of A, B, C, D,--- can be expressed in terms of the principal 
minors follows immediately from the fact that the elements of C, D,--- 
appear linearly and from the fact that we have been able previously to 
express the elements of A and B in terms of principal minors. 


Let us denote by J. the sum of all the principal minors of order « of 
the determinant A of order ». That the » sums thus obtained are 


independent becomes evident if we put a; — 0 (¢ + /) and observe that 
the sums J, (@ = 1, 2,---,m) are then simply the n elementary 
symmetric functions of aj (¢ = 1, 2,---, m). 


In the previous section we obtained new determinants of the nth order 
by combining the determinants A, B, C, D,---. We shall now study 
the sums Tupy..- of the principal minors of these new determinants, where 
Tuy... denotes the sum of all such principal minors containing @ columns 
from A, 8 columns from B, y columns from C, and so forth. 

We shall first find a superior limit to the number of independent sums 
Leggy... Obtained by the combination of two or more determinants. For 
this purpose let us observe that J,.3,... may be obtained from A by successive 
applications of two different types of differential operators. The first type 
is of the form 

If the kth successive application is denoted by k! Dk, we have 
DaA=Inr, = In. 


The second type is of the form 


De = > Dw— 


i=1j=1 0 aij 
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If the 7th successive application is denoted by /! Da», we have 
2 l 
Dav In 15 Dav In-x = In—k-2, 2; Dav In-x In-k-1,1. 


The operator 


which is of the second type, applied successively gives 


Thus the sum J,,... may be obtained from A by first operating with D, 
successively n—(a+A8-+y7-+---) times, by then operating with Da 
successively 8-+y-+--- times, by then operating with Dy, successively 
y-+.--- times, and so forth. 

We shall next show that all Z,,,... are solutions of the system of differ- 
ential equations 


n 


of 


k=1 O ks 0 Urk 


= 1,2,---,%). 


The determinant J, = A is evidently a solution of every equation of (17). 
It is easy to see that the relation 


=0 
0 Ars 0 Ars! 


Rrs (Da A) — Da(Rrs A) = 
is true, since terms containing the second derivatives vanish. But R,s(A)=0, 
whence it follows that R,s(Dg A) = Rrs(In-1) = 0. Therefore, J,—-1 is 
a solution of the system (17) and by induction it follows at once that all 
I, are solutions of (17). 

Let us repeat the above process with the operator Dg, instead of D,. 
We have 


= 
Ark 


Since Rys (Je) = 0, we have Rys(Dav Ie) = 0. It follows that is 
a solution of (17) and by induction it is again easily shown that all J,s 
are solutions. 


9 
{=1 j=1 0 bij 
0 Dks byx 
ol al 
Rys ( Dav Tu) Dav (Rys Tu) Zz (der Dsk «| 
k=1L\ O ks Ark 
0 Te 
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A continuation of the above process shows that all J,3,... are solutions 
of (17). 

The one relation >” , Rrr(f) = 0 between the equations (17) is obvious. 
That the remaining »?—1 equations are independent for two or more sets 
of elements aj, by, cyj,--- is easily seen by writing out the matrix of the 
equation and putting in this matrix. 

The system (17) thus involves kn* variables, where i is the number of 
the determinants A, B, C,---, and contains n*—1 independent equations for 
k>1. It follows that for k>1 there cannot be more than kn*— (n*—1) 
independent sums J,3,.... 

It is interesting to observe that for k = 2 and n> 2 the number 
}(n*?+3n) of sums J,,, is less than the superior limit n*+-1 fixed by the 
system of differential equations. As a matter of fact, there are }(n°—3n-+2) 
other independent solutions of equations (17) in this case.* 

We shall now prove 

THEOREM VII. The sums 


(18) (a, = 0,1, 2, ---, 2; y= 0, 1,2, 3) 


form a complete and independent system for k = 3. 

It is evident that there are » sums of the type Je and 4x(n+1) of 
the type (8 == 1,2,---,m). Consequently, there are n(n+1) of the 
type 4n(n+1)—~n of the type }n(m+1)—(2n—1) of 
the type J,,3,. These minors make a total of 2n?+1. 

Let us now put aj = 0 (7 +) and investigate the resulting functional 
matrix of the sums (18) with respect to au, bij, cj. There exists in this 
matrix a determinant of maximum order which consists of the product of 
the determinants 


and of a determinant F' consisting of the functional determinant of all the 
remaining sums of (18) with respect to all the elements bj, oj (4 + 7) except 
bisa (4 = 1,2,---,nm—1). The three determinants (19) are different from 
zero since each is equal to the functional determinant of the » symmetric 
functions Of Ann- 

We shall prove that F is different from zero by showing that the 
coefficient of the highest power of a,» which occurs in F is different from 


*Cf. Stouffer, Proceedings of the London Mathematical Society, ser. 2, 
vol. 15, p. 222. 
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zero. This coefficient is the product of the determinant which corresponds 
to F when the order of A is n—1 and the functional determinant 
Loe, Tour, Loo2, Tos, Loe, Toos, Lon, » Lo,n—2,25 Lo,n—8,3] 


(20) 


0 n—1, (n,n—1 Cn—1, Dn,n—2; Cn, n—2 bn ba; Cn1; Din, Cin] 


After forming the partial derivatives in (20) let us put all bx&(i +k) 
and cx (i+k) equal to zero except Dj i41(¢ = 1, 2,---,n—1), 
1,2, ---, n—383), and ci414(¢ = 1, 2,---, m—1). A partial 
derivative of Ig, with respect to bnr, brn, Cnr OF Crn 18 then zero unless 
at least one of the determinants in J), contains elements from the 7th row 
and column and also from all rows and columns of higher number. It 
follows that (20) becomes the product of the determinants 


(21) 
[ Dur, Cnr Cra ] 
and the determinant 
(22) 


Loc, Jon, Looe | 


0 [ On, Cn,n—15 Ca—1, 2) 


Each of the determinants (21) and (22) is different from zero since, under 
the above conditions upon by and cj, each is the product of the elements 
of its principal diagonal, all of which are different from zero. Thus F' does 
not vanish when the order of A is n, provided it does not vanish when 
the order of Ais m—1. But it is easily verified that F’ is different from 
zero for n= 3. Hence, by induction the sums (18) are proved to be 
independent. Since the set (18) contains the maximum number of independent 
sums, it is a complete system. 

The nature of the determinants (19) and the fact that (20) reduces to the 
product of the elements of its principal diagonal prove the independence of 
the sums (18) and 


for k = 4, 5,---, and any v, provided that they are independent for n = 3. 
This latter fact is easily verified. 

Since the set (23) contains x? sums for 06 = 1, n* sums for ¢e = 1, and 
so forth, we have proved 

THEOREM VIII. The sums 


(24) (y = 0,1, 2,3; d+e+ 


form a complete and independent system for any number k of combined 
determinants of any order n. 


= 1,2, ---, n—2) 

= 1) 
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IV 

Dickson* has determined the types of general homogeneous polynomials 
which are expressible as determinants with linear elements. We shall now 
show that his main results follow very easily from our knowledge of the 
independence of sums of principal minors of determinants. 

Let f(2,, +++, Ze) represent a general homogeneous polynomial of 
degree n in k variables. We assume that the coefficient of the nth power 
of one of the variables, say 2,, is unity. There is no loss of generality 
here since the coefficient can always be made unity by a linear trans- 
formation. If it is possible to express f as a determinant with linear 
elements, a combination of rows and columns will always bring that deter- 
minant to the form 


| 
Az, + des + Aon tet bena'st - 


j ! | 


The expansion of (25) gives a polynomial in which each coefficient 
after the first is one of the sums J,,,... and, conversely, each sum 
(a, 8, v,--- =0,1,2,---,m) appears as a coefficient, the number 
of subscripts in J,,,,... being /—1. Since this polynomial must be exactly /, 
we have a proof of 

THEOREM IX. A general polynomial of degree n in k variables cannot be 
expressed as a determinant with linear elements if it has in addition to 2’ 
more than (k—2)n*?+1 terms if k>3 or more than 4(n?+3n) terms 
if k = 3. 

This is Dickson’s first theorem except that he included the case k = 
in the general theorem. For k = 3 the limit is }(n?+ 3m) which for n > 
is smaller than n?+-1, the limit for k = 3 in the general theorem. 

THEOREM X. Any determinant of order n whose elements are linear homo- 
geneous functions of x71, X2,-++, ® can be expressed in a canonical form 
involving not more than (k—2)n?+2 parameters for k > 3 and not more 
than 4(n®+3n)+1 parameters for k = 3. 

This is Dickson’s generalization of his first theorem except for the 
case k = 3. Its proof follows immediately from the maximum number of 
independent sums which appear as the coefficients in the expansion of the 
determinant. The additional parameter arises from the fact that the 


*These Transactions, vol. 22 (1921), pp. 167-179. 
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coefficient of 2” was assumed to be unity in the expansion of (25). If the 
expansion of the determinant does not contain the nth power of one of the 
variables, a linear transformation, independent of the coefficients of the form, 
will introduce such an nth power. The above proof then applies. The 
inverse of the linear transformation will restore the determinant to its original 
form without increasing the number of parameters. 

Our knowledge concerning the particular sums J,,3,,. .. which are independent 
gives immediately a limitation as to the forms which we may hope to express 
as determinants with linear elements. For n = 2 or 3, the coefficients J,,5,,.. 
are not all independent for k > 4. For x > 3 they are not all independent 
for k > 3. Thus, we have 

THEOREM XI. Unless k = 2 or 3, n any integer, or k = 4, n = 2, 3, the 
general form is not expressible in determinantal form. 

This is Dickson’s second theorem. 

Within the limitations just stated there is a one-to-one correspondence 
between the independent coefficients J,,,... of the expanded form of (25) 
and the coefficients of the general form. Consequently, we have 

THEOREM XII. Jn the field of all compler numbers the general binary 
form, the general ternary form, the general quaternary quadratic form, and 
the general quaternary cubic form can be expressed in determinantal form. 

This is Dickson’s fifth theorem, except that Dickson shows that the 
theorem is true for every binary form, every ternary form, and every 
quaternary quadratic form, regardless of the generality of the form. 

The fact that the general quaternary cubic form can be expressed in 
determinantal form with linear elements leads at once to a theorem proved 
by H. Schréter* and also by Cremona;t 

THEOREM XIII. A sufficiently general cubic surface can be generated by 
three projective bundles of planes. 

For the determinantal form of the cubic equation says at once that 
there exist three parameters 4, ~,v such that the three equations 


Ala = 0 (¢ = 1, 2, 3) 


are satisfied by the codrdinates of all points of the surface, where /;j are 
expressions linear in %2, 2%. 

* Journal fiir Mathematik, vol. 62 (1863), p. 265. 

+ Ibid., vol. 68 (1868), p. 79. 
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A NECESSARY AND SUFFICIENT CONDITION 
THAT TWO SURFACES BE APPLICABLE* 


BY 


W. C. GRAUSTEIN anp B. 0. KOOPMAN 


It is well known that, in order that two surfaces be applicable, it is 
necessary that a map of the one upon the other exist so that geodesics 
correspond to geodesics and total curvature be preserved. It is also 
a familiar fact that neither of these conditions is alone sufficient. The 
primary purpose of this paper is to show that the two conditions taken 
together are sufficient, i. e., to prove the theorem 

If two surfaces can be mapped geodesically so that total curvature is pre- 
served, the surfaces are applicable. 

1. The point of departure for the proof is Dini’s theorem to the effect 
that, if two surfaces correspond by a geodesic map, then (a) each is mapped 
isometrically on the other or on a surface homothetic to the other, or 
(b) the two surfaces are surfaces of Liouville, whose linear elements can 
be put simultaneously into the forms 

S. 


1° 


= (U+V) (du*?+ dv’), 
(1) 


du? 


where U and V depend, respectively, on uw and v alone, and corresponding 
points have the same curvilinear coérdinates. 

When we demand, further, that the geodesic map preserve total curvature, 
the surfaces in case (a) are obviously applicable. Case (b) is disposed of 
by the following lemma: 

If two surfaces of Liouville with linear elements of the forms (1) have the 
same total curvature in corresponding points, they are surfaces of constant 
curvature. 

For, it follows then that the surfaces are applicable, though not, it is 
to be noted, by the correspondence established by equations (1). 

To prove the lemma, we compute the total curvatures K, and Ky, of 
S, and S, by means of the Gauss formula 


* Presented to the Society, December 27, 1922. 
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1 
2(U-+V)8 


1 


K, v"—(U+V) (U"+V")], 


K, = [—2 (U*V*+ UV*)U" + 


+2(U*V?+ 
Setting K, = Kz, we have 
2(U+V—U? U" + v3—2)U"”" 
(3) 
U3 +2)V" = 0. 


By means of the substitutions 


(4) 


(3) becomes 


,, dX 
(at + (3xy?+ y®—2) X 
ar 


ly 


—(32%y+ 2) = 


where X and Y depend, respectively, on x and y alone. 
Differentiating (5) four times with respect to x, we get 


+ (2— 3y>— 5ry*) £ 


dat = 0. 


(a 
Since x and y are independent variables, it follows that d*X/dx* = 0. 
Similarly, d*Y/dy* = 0. Hence 
X = azz’, Y= 


Substituting these values of X, Y in (5) and equating collected coefficients 
of «” y" to zero in the result, we get 


Hence, by virtue of (4), 


= a,U+ a, 


(6) 
= a, +a43(V°+1). 
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(2) 
y=V, X=0% yY=r", 
(5) = 

dy 
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On substitution of these values in (2), we find that Kj = — 74s, so that 
K, and Ky, are constant, and our proof is complete. 
2, Equations (6) form a special case of the equations 


a,U+ 


= a, V—a,V*+ a3 


For these more general values of U and V, 


Conversely, if the curvature of either of the surfaces S,, S; is constant, 
U and V must be defined by equations of the form (7). In proving this, 
there is no loss of generality in assuming S, to be the surface of constant 
curvature, for the relationship between S, and S, is reciprocal. Accordingly, 
we set Ky = — ; a3 in (2), obtaining the equation 


(lg 


(U+ V) Vv") — (Ur + Vv’) (0+ == 


On application of the substitutions (4), this reduces to 


aX dy 
(8) 2) ag (x+y)’. 


Differentiating twice with respect to x, we get 


= 6a3, whence = ax’. 
Similarly, 


and Y= &+hyth 


Determining the coefficients in X, Y by substituting in (8), we come out 
with equations for U and V of the desired form (7). Incidentally we have 
also proved that the only surfaces which can be mapped geodesically 
on a surface of constant curvature are surfaces of constant curvature— 
Beltrami’s theorem in a generalised form. 


(7) 
a a 
4 4° 
8 
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If in (1) U and V are replaced by U'-+h and V—h, where h is an 
arbitrary constant, S; is unchanged, but S. is replaced by a one-parameter 
family of surfaces. The same substitutions in (7) leave as unchanged and 
replace ad) by agh*+ Consequently, we have 


as (lo + ay h + (lg Ag 
Ks 
4 


Thus corresponding to a given surface S,, that is, for a given set of values 
for the a’s, there exist three surfaces S., in general distinct, of the same 
constant curvature as S,;, namely those corresponding to the three roots 
of the equation 


as as h?+- + Ay+ a, = O. 


An exception arises in case S, is a developable (a; 0); there then exist 
among the surfaces S, at most two developables. 
HARVARD UNIVERSITY, 
CAMBRIDGE, Mass. 
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EXTENSIONS OF RELATIVE TENSORS” 


BY 


OSWALD VEBLEN anp TRACY YERKES THOMAS+ 


1. Introduction. This paper is intended as an addendum to our 
previous paper on the Geometry of paths.t It contains the deduction of the 
formulas for the covariant derivatives and the higher extensions of relative 
tensors. These formulas are calculated by a process entirely analogous 
to that by which the corresponding formulas for ordinary tensors are 
obtained. Although the formulas are rather obvious the only one which 
we have observed in the literature is that for the covariant derivative of 
a relative tensor of the second order which was used by Einstein.§ 

2. Relative tensors. Let us denote by 


(2.1) = 


an arbitrary analytic transformation of the coérdinates (2’,-- 


can be written in the inverse form 


-, 2”) which 


(a',---, a”). 


Since the inverse transformation (2.2) exists, the jacobian of the trans- 
formation (2.1) 


| | 


Ox | 


(2.3) A 


does not vanish identically. As we shall later require the derivative 0 A/dzx* 
of the jacobian A we here note that it is given by 


OA 


(2. 
4) 


* Presented to the Society, April 19, 1924. 
7 National Research Fellow in Mathematical Physics, University of Chicago. 
} These TRANSACTIONS, vol. 25 (1923), p. 551. 
§ A.Einstein, Zur allgemeinen Relativitatstheorie, Sitzungsberichte der PreuBischen 
Akademie der Wissenschaften, 1923, p. 32. 
373 


(2.2) 
axe 
ox" ax” Az 


OSWALD VEBLEN AND T. Y. THOMAS [July 


A set of functions 7)j...~° will be said to form a relative tensor of 
weight A if it transforms according to the equations 


K 
A 


when the variables are transformed by equations (2.1). If K = 0, 77"; 

is an ordinary tensor. In case K Ti" is said to be a tensor 
density and is then usually written see The justification of the name 
tensor density lies in the fact that the law of transformation of the integral 


colm---n 


extended over a definite »-dimensional region approaches more and more 
closely the law of transformation of an ordinary tensor Thx as the 
region of integration closes down on a point. Thus the tensor density /"";;" 
represents a sort of tensor 7"";;" of weight zero per unit of codrdinate 
volume dV. 

We easily see that the algebraic processes of addition, multiplication, 
and contraction of ordinary tensors will also hold for the case of relative 
tensors. 

3. Covariant differentiation or first extension. Let the relative 
tensor 7'/"";." of weight K which is referred to the codrdinates (', ---, 2”) 
be denoted by 4/"";" when referred to the system of normal codrdinates 
(y',---, y") which are determined by the coérdinates (x',---, 2") and a point 
(q',---,q"). We shall show that 


\ 


(3.1) Tip % +k 


0 
Oy? 


where the derivative is evaluated at the origin of normal coérdinates, 
defines a set of functions Ts" of (z',---,a") which are the components 
of a relative tensor of weight K. The relative tensor 77x.) will be called 
the covariant derivative or first extension of the relative tensor Tix 
We denote by (y',---,y”) the system of normal coédrdinates determined 


by the codrdinates (a',---,2”) and the point (q',---,q"). Then 


29 i i 
(3.2 y ay 


on 

374 
“ls m S o,f 3 4 
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n 
where the a’s are constants. Furthermore denoting by ¢;;-:-;" the relative 


ylm-+++n 


tensor 7ij...«. When referred to the normal codrdinates (y',---, y”) we have 


Im-+-n K O ay” dy" 
dy" oy ay 


In view of (3.2) the derivatives in (3.3) are constants and hence 


dy? dy" dy” ay? dy’ dy! 


Evaluating at the origin of normal codrdinates we obtain 


lm- 


(3.5) Ty. 


K 
A 


batt Oa” 


which shows that the set of functions 77/"";:", constitute a relative tensor 
of weight A contravariant in the indices (/,m,---,) and covariant in 
the indices (7, 7,---,k,p). 
We note that formulas analogous to those for the sum and product of 
two ordinary tensors likewise apply for the case of relative tensors. 
To obtain the explicit Pe involving the 7’s and their derivatives 
n 


for the relative tensor .-k,p We make use of the equations 


K 
Og Or Oy 0 yf 


Differentiating these equations and evaluating at the origin of normal 
coérdinates we obtain 


a 
Ty. Vp ++ +++ Ti. 


Ti 
Oa? 


7 
— Taj... lip - Ti Vip — K 


This is the general formula of covariant differentiation. By using the 


tensors x. and Egy... (ef. The geometry of paths, loc. cit.) the 
form (3.7) may be abbreviated. Thus 


Oa? 


im - atp 
‘yy 


ot 
(3.4) 
dy?” 
ax! dum da? Oa” 
» 
(0.8) n r 
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e 


4, Higher extensions. By a process similar to that employed for the 
covariant derivative we may show that 


0 


defines a set of functions of (2',---, 2") which constitutes 
a relative tensor of weight K. The relative tensor 7'/"\i’»...9 will be 
spoken of as the rth extension of the relative tensor Ta weovited that 
there are r indices (p,---,q). In case r—1 the extension reverts to the 
covariant derivative which we have already considered. 

From its definition by means of the equations (4.1) we see that the 
extension ‘k,p---q symmetric in the indices (p,---,q). Thus 


im: lm---n 
(4.2) Ty. = Ti. u- 


where (w,---,v) denotes any permutation of the indices (p,---,q). 

The formulas for the extension of the sum and product of two relative 
tensors are similar to the corresponding formulas for the extension of the 
sum and product of two ordinary tensors. 

General formulas of extension (7 > 1) may be calculated by the same 
process as that employed in the calculation of the general formula of 
covariant differentiation (3.7). The formula for the rth extension 

lm.-- 
of a relative tensor 7)j".;;" of weight K +0 involves the formulas for 
the first x extensions of 7jj"\." considered as a tensor of weight zero. 
For we have 
‘ 
where 

Oy OY os 


Differentiating (4.3) and evaluating at the origin of normal codrdinates 
lm. - lm 
Im- 
(4.4) + S (Apg ) +-- 


K 
+8 (Apgr --u Kk v -uv Lij---k 
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where 


K 


and S( ) denotes the sum of the different terms obtainable from the 
one in parenthesis by forming arbitrary combinations of the subscripts 
The expressions are given by the ordinary 
formulas for the rth extension of a tensor 7y".:" of weight zero, where r 
is the number of indices (p,---,q), but these expressions do not in general 


constitute a tensor. The quantities £. y have the values 


ap? 


a a 


KS 
re apqr KS ( ye) ( nie) 


ap 


in which Ij ete. are definite functions of the 7%, and their derivatives 
The geometry of paths, loc. cit., p. 561). Thus we may write 


ly 
(4.5) Ti = KT 


in place of the formula (3.7). 

Any formula of extension of a relative tensor Tj x of w eight K may 
be obtained by substituting the proper values of and Ax .. 
in (4.4). A few particular expressions 7//”""%.... are given in the Senate 
of paths as formulas of extension of tensors of weight zero and these when 
substituted into the equations (4.4) together with the above values of i. 
give complete formulas in terms of the Tix and their derivatives for the 
extensions of certain relative tensors of weight K. 
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GEOMETRIES OF PATHS FOR WHICH THE EQUATIONS 
OF THE PATHS ADMIT A QUADRATIC FIRST INTEGRAL* 


BY 


LUTHER PFAHLER EISENHART 


1. A geometry of paths for a general space of m dimensions is based 
upon the conception that the paths are fundamental entities of the space, 
They are the integral curves of a system of differential equations 


(1.1) 


where = 1, ---, m) are the codrdinates of a point, are 
functions of the z’s and s is a parameter peculiar to each path. An 
important class of these geometries, which admits the Riemannian geometry 
as a sub-class, is that for which the equations (1.1) admit a quadratic first 
integral 


(1.2) 


where g; (= g,;) are the components of a covariant tensor. Veblen and 
Thomast have considered the problem, given a set of /’s to determine 
whether equations (1.1) admit an integral (1.2); they have shown that its 
solution involves only algebraic processes. In this paper the converse 
problem is solved, namely to determine the /”’s so that (1.1) shall admit 
a given first integral (1.2); also the more general problem when the first 
integral is of the form 


ea dai 


= COnst., 
Iij ds ds 


where 9 is a vector and the integral is taken along the path. 
2, General formulas. If we puta! = thus introducing 
a new set of codrdinates, equations (1.1) become 


* Presented to the Society, March 1, 1924. 
+ These Transactions, vol. 25 (1923), pp. 599-608. 
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(2.1) 


where 


Suppose now that Gj are the components of a symmetric covariant tensor 
of the second order, that is 


0 at 
(2.3) Jap = 95 
Ox 


then by means of (2.2) it can be shown that the functions g,,,, defined by 


a a 
(2.4) Vie Pin Vins 
are the components of a covariant tensor of the third order. As thus 
defined g,,,, is a generalization of the first covariant derivative of 9, the 
ordinary covariant derivative for a Riemann space being given by (2.4) 


the Christoffel symbol of the second kind 


formed with respect to the fundamental form of the space,* and similarly 
for If. 

The components g/ of the contravariant tensor associate to g,; ave given by 


when 7° is replaced by ; 


(2.5) Ie = 


j 
where 


(2.6) dj = 0or 1, ori=j. 


As thus defined, g¥ is the cofactor of g,; in the determinant g = |g,,| 
divided by g. 
If we put 
09; OG 09g;; 
2 \dx/ + 0 x* 
and 


(2.8) = g§*[ij, «], 


Jk 


* Cf. Proceedings of the National Academy of Sciences, vol. 8 (1922), p. 21. 
27 


379 
d? da’ La’? 
Kk a Jj yk 
1B hed ap ad 
Ox Ox Ox 
e 
n 
y 
d 
e 
is 
st 
1g 


380 L. P. EISENHART [July 


then [ij,k] and at are the Christoffel symbols of the first and second 


kinds respectively formed with respect to the g’s. 
If we write also 


(2.9) = gal, Te = 
then from (2.4) and similar equations we have 
(2.10) + — = KI — 
From (2.3) we have by differentiation and suitable operations* 


(2.11) 


where he is formed with respect to the g’’s. Subtracting this equation 


from (2.2), we obtain 


Hence if we put 
(2.12) 
the functions ¢c,,* are the components of a tensor contravariant of the first 


order and covariant (and symmetric) of the second order. In consequence 
of (2.8) and (2.9) we have from (2.12) 


(2.13) = Liz, + 
Hence (2.10) becomes 


(2.14) j + Dix, “faa — 


If we add to this equation the two equations obtained from it by per- 
muting 7,7, k cyclically, we obtain 


(2.15) + Iix,i + j = —2 (cijx + Cjki + Ckij ) 


* Cf. Bianchi, Lezioni, vol. 1, p. 64. 
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3. Quadratic first integrals. When we express the condition that (1.2) 
be a first integral of (1.1), we obtain 


dé do at 
Iij,k ds ds ads 


Since this must be satisfied identically, we must have 


(3.1) + = (i,j,k=1,---,n). 


The consistency of these equations is the necessary and sufficient condition 
that (1.2) be a first integral.* 
We consider also the case when (1.1) admits a first integral of the form 


(3.2) const., 


where % are the components of a vector and the integral Jv dx* is taken 
along the path in question. Proceeding as above, we obtain the equations 


Conversely, if these equations are consistent and yield a tensor g, and 
a vector Qa, then (3.2) is a first integral of the corresponding equations (1.1). 
It can be shown7 that if we put 


(3.4) = 4 vi) 


and define a parameter s for each path by means of the equation 


(3.5) 


* Veblen and Thomas, loc. cit., pp. 599-608, give a complete treatment of the question 
of the consistency of equations (3.1). 
+Annals of Mathematics, ser. 2, vol. 24 (1923), p. 376. 
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the equations 


(3.6) 


ds* Yds ds 


define the same curves as (1.1). 
If we denote by Ii,% the covariant derivative of g,, given by the equation 
obtained by replacing is by Tj in (2.4), we have 


Now (3.3) becomes 
(3.8) + + j = 0, 


and (3.2) reduces to the form (1.2) by means of (3.5). 

The change (3.4) in the /’s means a change in the affine connection of 
the space, but not in the paths themselves. Hence we have the theorem 

When the equations of the paths of a space admit a first integral of the 
Sorm (3.2), by a change in the affine connection but not in the paths them- 
selves the equations of the paths can be given a form which admits a first 
integral of the form (1.2). 

When, in particular, 


(3.9) = ~ Gg Pre 


equations (3.3) are satisfied. This is Weyl’s geometry, for it follows 
from (2.10) that 


Consequently, if the paths are taken as fundamental rather than the affine 
connection, it follows from (3.4) that if in (1.1) we take 


(3.11) rk — + (dF gi + gi 
we have a geometry of the space with the same paths, whose equations 
admit the first integral (1.2). 

Returning to the consideration of (3.1), we observe that from this equation 
and (2.15) we have the theorem 


* Cf. Space, Time and Matter, p. 125. 
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A necessary and sufficient condition that equations (1.1) admit a quadratic 
first integral (1.2) is that the tensor c;,, defined by (2.13) satisfy the conditions 


(3.12) Cun + Cini = 0. 


4, Determination of geometries of paths with quadratic first 
integrals. Suppose that we have any symmetric tensor g,; anda covariant 
tensor of the third order, a,,, symmetric in ¢ and j. If we define a tensor, ¢;,,., 
by means of the equations 
(4.1) 2 — — (i,j, = 1,---, 0), 
it is by definition symmetric in 7 and j, and satisfies (3.12). If we define 
a set of ’s by means of (2.13) in which [7j, k] are formed with respect 
to the given tensor g,;, and ¢,, is given by (4.1), it follows from the theorem 
at the end of § 3 that the equations (1.1) admit the corresponding first 
integral (1.2). 

Consider, conversely, the case when a geometry of this kind is given. 
Then the ¢,, as given by (2.13) satisfy (3.12). This condition is met, if we 
determine the components a,,, of a tensor such that (4.1) hold for ¢,,, known. 

For each set of values of 2, 7, k all different, there are two equations (4.1), 
which are equivalent to 


Gini = dix +); 


where 


the second and third expressions being equivalent because of (3.12). 
From (4.3) follow the identities 


djix + dey = 0, 
(4.4) 
djix + + = 0 k= N; k ). 


In the general case, that is when c,, do not satisfy any conditions other 
than (3.12), there are, in consequence of (4.4), 2(2— 1) (2—2)/3 independent 
equations of the type (4.2). 


y 
n 

(4.2) 
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From (4.1) and (3.12) we have also 


where 


(4.6) 


In the general case there are n(m—1) independent equations of this type. 

From (3.1) we have g,,; = 0, so that from (2.10) we have I,,, = [77, 7] 
and consequently ¢,;, = 0. In this case (4.1) vanishes identically and 
consequently the components a;,, are not determined. Hence there are 
n(n—1)(n+1)/3 independent equations for the determination of the 
n*?(n-+1)/2 components of Oy Consequently there aren (n-+1)(n-+ 2)/6 
components arbitrary: they are one for each case where 7, 7, k are different; 
one where two are different; and all of the type a,,. Hence we have the 
theorem 

A tensor g,, and a tensor a;,,, symmetric ini and j, determine a geometry 
of paths whose equations admit the corresponding first integral (1.2); con- 
versely, if a geometry is given whose equations admit a first integral (1.2), 
n(n+1)(n+2)/6 of the components Aj, are arbitrary and the others can 
be found directly. 

PRINCETON UNIVERSITY, 
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A GENERAL MEAN-VALUE THEOREM* 


BY 


D. V. WIDDER 


In a paper published in 1906+, Professor G. D. Birkhoff treated the mean- 
value and remainder theorems belonging to polynomial interpolation, in which 
the linear differential operator u™ played a particular rdle. It is natural 
to expect that a generalization of many of the ideas of that paper may 
apply to the general linear differential operator of order m, and the author 
is attempting such a program. This generalization throws fundamentally 
new light on the theory of trigonometric interpolation. 

A very elegant paper by G. Pélyat has just appeared treating mean-value 
theorems for the general operator in a restricted interval. It is the special 
aim of the present paper to develop a general mean-value theorem, and 
to show how it can be specialized to obtain Pélya’s results. 

We consider a linear differential expression of order n, 


Lu = (x2) +h (2) 


where J,(x), are continuous functions, and w(x) is con- 
tinuous with its first (1—1) derivatives, the nth derivative being piece- 
wise continuous. All functions concerned are real. It is the purpose of 
this paper to obtain a necessary and sufficient condition for the change of 
sign of Zw in an interval in which w vanishes (m-+1) times. 

More generally, the (x -++1) conditions implied in the vanishing of wu may 
be replaced by an equal number of conditions involving also the derivatives 
of wu. Let x, 71,-°--, 2n be points of the closed interval (a, 6), which 
points need not be all distinct, and let ko, ki,---, kn be zero or positive 
integers not greater than n—1. We take then as »+ 1 conditions on u the 
relations 


(A) uk) = 


* Presented to the Society, May 3, 1924. 

+ G. D. Birkhoff, General mean-value and remainder theorems, these Transactions, 
vol. 7 (1906), pp. 107-136. See also Bulletin of the American Mathematical 
Society, vol. 28 (1922), p. 5. 

tG. Pélya, On the mean-value theorem corresponding to a given linear homogeneous 
differential equation, these Transactions, vol. 24 (1922), pp. 312-324. 
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Here u‘*)(a) denotes the kth derivative of u, and w(x) is the same as 
u(x). We assume that no two of these equations are identical. 

Now let 7%, w2,--+, Un be m linearly independent solutions of the homo- 
geneous equation 


(1) Lu = 0. 


For definiteness take them as the principal solutions for the point a; that 
is, solutions satisfying the conditions 


m(a)=0, wla)=0, ---, = u*—Y(a) = 1, 


uz(a) = 0, uz(a) = 0, (a) = 1 (a) = 0, 
= 1, u,(a) = 0, (a) = 0. 
We consider also the non-homogeneous equation 
Lu = g(z), 
where g(x) is piece-wise continuous in (a, b). The general solution of (2) 
is now obtained by Cauchy’s method. Determine a solution of (1) which 
together with its first »—2 derivatives vanishes at a point ¢ of (a, b), while 


the (n—1)th derivative has the value unity at that point. Denote the 
function by g(x, t). It satisfies the n conditions 


(3) g(t, t) =0, g(t, t) =90, g”*(t, t) = 0, (4, t) =1. 


Here the differentiation is with respect to the first argument. 
It is known that the general solution of (2) may be written in the form 


b 
u(r) = +eu+--- +9(x,t) p(t)dt. 


Here 1, ¢2,--+, (mn are arbitrary constants, and the sign before g(z, ¢) is 
to be taken positive if ¢<a and negative if ¢>«x.* Considered as a 
function of ¢, g(x, ¢) satisfies the equation adjoint to (1), 


d 
dx 


d”y 


—1 
= +1?" v)+ +h,v0 =0. 


*See for example D. A. Westfall, Dissertation, p. 16. 
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If now we express the fact that «(2°) satisfies the n+ 1 conditions (4), 
we obtain n-+1 equations 


a 


(¢ = 6,1,3,---, 2%). 


Eliminating ¢,, c2,---, Cn from these equations, there results an equation of 
the form 


b 
(5) == ©, 


where A(t) is the determinant 
| + 4g (x,, t) uy (x;) us? (x;)++- (x;)| 


Denote the cofactors of the elements of the first column by Ao, 4:,---, An, 
so that A(t) takes the form 


A(t) = (ai, t) di. 
4=0 


It is evident that A(t) depends in no way on the choice of the linearly 
independent solutions u,, we, ---, un, but merely on the position of the 
points 2%, 

Now let us suppose that 4(¢) is not identically zero. Then the function wu 
satisfying the conditions (A) can not be a solution of (1) unless it is identi- 
cally zero; for a necessary and sufficient condition that there exist a solution 
of (1) not identically zero and satisfying the conditions (A) is precisely that 


Then g(t) is not identically zero, and we have at once from (5) a sufficient 
condition that Lu change sign in the interval (a, b), namely that 4(¢) should 
be a function of one sign in that interval. We may in particular take 
a and b as the two points of the set x, z,,---, 2% Which are farthest apart, 
and thus be assured that the change of sign of Lu occurs between these 
two points. 

The condition is also necessary. Suppose that any function « with the 
required degree of continuity which is not identically zero and which satis- 


uly 
at 
), 
), 
(4 = 0, 1,---,%). 
e 
1 
(¢ = 0,1,---,%). 
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fies the conditions (A) is such that Lw changes sign between the extreme 
points. It is desired to show that A(t) is a function of one sign not identi- 
cally zero. 

We note first that the 4; are not all zero; for if they were, there would 
be a solution wu of (1) satisfying the conditions (4). This is impossible 
since Lu must change sign by hypothesis. In order to prove that 4(¢) does 
not vanish identically we must investigate its structure more closely. 

If we denote by 7, vs, ---, Un the solutions adjoint to u,, we,---, Un, Wwe 
may write g(x, t) as follows:* 

n 


g(a, t) = 2 ti (x) v(t). 


A(t) then takes the form 


If 4(t) were identically zero, each summation in the above expression would 
be zero, since v,, v2,---, Un are linearly independent. Now by taking ¢ in 
all possible positions with respect to the points a9, 2, ---, Zn, Various 
combinations of signs in each summation are obtained. 
It would follow then that 
J 


(x,) 4, = 0 m3 = 1, 2,---, 2). 
Not all the 4; are zero. Suppose that 4, +0. Then it would follow that 
u(x.) = 0 


The Wronskian of %, ws, ---, u, Would vanish at the point 7», contrary 
to the assumption that u,, w2,---,%» are linearly independent. A(t) can 
not therefore vanish identically. 

Suppose now that 4(¢) changes sign between the extreme points. It is 
then possible to choose a continuoust function of one sign y(t) such that 


b 
(6) g(t) A(t)dt = 0. 


* See for example, Darboux, Théorie des Surfaces, vol. 2, p. 106. 
+ Indeed @ may be a piece-wise continuous function made up of straight lines parallel 
to the x-axis. 
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Now consider the equations (4) in which g(t) is replaced by g(¢). From 
these n-+1 equations, pick out that set of » equations which has 4, for 
its determinant. Since 4,, +0, the set of equations has a unique solution 
in (1, C2,***, Cae If we substitute this solution in the equation 


*b 
(x) = +entn(w) +4 + t) p(t) dt, 
e/a 


we obtain an equation of the form 


(2) = (x, t) p(t) dt, 
e/a 

where G(a,t) may be identified with the Green’s function corresponding to 
the boundary conditions (A).* 

It is seen that 

G“” (am, t) = A(t). 
It follows from (6) that 
u™ (am) = 0. 


By its very definition wu (a) is seen to satisfy the remainder of the con- 
ditions (A). Lu must therefore change sign. But Lu is g, a function 
which does not change sign. We have thus completed the proof of the 
following 

THEOREM I, A necessary and sufficient condition that Lu change sign 
in an interval in which u (having the required degree of continuityt and 
not identically zero) satisfies the conditions (A) is that A(t) be a function 
of one sign not identically zero in that interval. 

G. Pélya obtains certain theorems concerning the vanishing of Lu. We 
may obtain these results from Theorem I. 

With Pélya we say that the property W holds for the operator Lu in 
an open interval (a,b) if there exist solutions of (1), J, he,---,n-1, such 
that the following functions do not vanish in (a, b): 


hy he 


Wi hy, We = W (ha, he) hi 
1 


| hy he 


*C. E. Wilder, these Transactions, vol. 18 (1917), p. 416. 
7 The restrictions on the nth derivative of u might be made lighter as is done in 
Birkhoff’s paper, loc. cit. 


ly 
1e 
ld 
le 
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Pélya considers the special case of the conditions (A) which involves the 
vanishing of «(a) at points of the interval that are distinct or coincident, 
Consider ry points 


Ly +++ ay, 
Suppose that u(a) vanishes m; times at a point a; (m<—n—1): 


u (xi) = (xj) = --- = (a4) = 0, t—1,2,---,r; 
(B) 


> m=nt+i. 


i=1 


If 21, 22,---+, 2, lie in the interval in which the property W holds, A(t) is 
a function of one sign. To prove this we investigate the structure of A(t). 
In each of the r —1 intervals, A(¢) is a solution of the adjoint equation, 
continuous with its first m derivatives. It is only at a point 2; that 
a discontinuity may occur. Such a discontinuity is caused by a change 
in the ambiguous sign before g*® ( xi, t) as t passes over 2;. However, it 
is only when g“* (2;, x;) is not zero that such a discontinuity is introduced. 
We can now show that at a point 2; where u(x) vanishes m; times A(#) 
is continuous with its first » —m;—1 derivatives. 

It is known* that 


g(x,t)| =D 
it 


da" at” 


2 t=2 


At xj, # may equal m;—1, so that 


=a” 
It follows then that A(¢) is continuous at 2; with its first n—m;—1 
derivatives. 

Now if t<a, i = 1, 2, ---, vr, then all the ambiguous signs in 4 (f) 
are positive; if t>a;, i = 1, 2, ---, r, they are all negative. In either 
case 4(t) is identically zero for the interval considered, since 


uj (24) A, = 
i=0 


* See Schlesinger, Lineare Differential-Gleichungen, vol. 1, p. 63. 
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the expression on the left being a determinant with two columns equal. Now 
since 4 (t) is continuous with its first » —m,—1 derivatives at z,, it follows 
that A(t) has n—m, zeros at z,. It has n—~m, zeros at zy. 

Suppose first that 


mM, = My = = 1, 


Then 4(¢t) is continuous throughout with its first »—2 derivatives and 
has »—1 zeros in each of the points 2, and 2,. Suppose it were not 
a function of one sign. It would have in all at least 2%—1 zeros in the 
closed interval (2, 2). 

Now recall that I/(v) can be “factored” as follows: 


(—1" d Wi W a 


= dz dx dc WaWes dz 
where Wo = 1.* Each of the quantities Wo, Wi, ---, Wn does not vanish 
in the interval (21, zn+1), so that we may apply Rolle’s theorem. If A (1?) 
vanished times in (2, the function 


y(t) = W.Wo dx WiWs dx dx A(t) 


would change sign at least » times inside the interval. But M(A(t)) is 
identically zero, so that w(t) is constant in any interval in which it is 
continuous. Hence w(t) can change sign only at the points x2, 73,---, Xn, 
where it is discontinuous. There are only »—1 such points, so that it 
must be concluded that 4(¢) is a function of one sign in (21, 2n+1). 

In order to obtain the proof in the general case we shall need the 
following 

Lemma. Jf a function f(a) is continuous in an interval in which f" (x) 
is continuous except for l finite jumps, and if f' (a) can have at most N zeros 
in the interval, f(a) can have at most N+1-+1 zeros there. 

Proof. If f(x) had N-+/-+ 2 zeros, f’(x) would have N+/+1 zeros 
and discontinuous changes of sign. At most / of these can be discontinuous 
changes of sign, and f’ (x) would have N+ 1 zeros contrary to hypothesis. 

Now denote by s the number of integers mz, ms, ---, mr—1, Which are 
equal to 7. Then 


= r—2. 


*See Schlesinger, loc. cit., p. 58. 
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Define a function 4''(¢) by the equation 


Now’apply the lemma to A” in each of the »—s,—1 intervals in which 
it is continuous. We may evidently treat all these intervals simultaneously 
and take / equal to the total number of discontinuities of _ , viz., r—2., 
N must be zero since A” ” is constant where it is continuous.* We con- 
clude that A”~” can vanish at most r—1 times. 

Now apply the lemma to Here N= r—1 and 1=r—s,—2. 
Then A” can have at most 2(r—1)—~s, zeros. A” can have at 
most 3(r—1)—2s,—s,. Proceeding in this way we see that A(t) 
has at most 


(n —1)(r —1) —(n— 2) 8, — — 3) — - 
zeros. Now A(t) has 


N—m+n—- Mm = 


zeros at the end points x, and 2,, and this number is precisely equal to 
the maximum number of zeros A(t) can have, since by virtue of (7) 


(n —1)(r—1) —(n— 2) 5, — (n — 8) --- = 


The proof is thus complete that A(t) can not change sign. 

Before proceeding to the converse of this theorem, let us draw several 
further inferences. 

THEOREM II. Jn an interval in which the property W holds, the coefficient 
of li (aj, t) in A(t) can not vanish, i= 1,2,---,r. 

For if it did A(t) would be continuous with its first 1 — m; derivatives 
at 2;, and by means of the lemma a contradiction would be reached as 
before. 

CoROLLARY. No solution of equation (1) can vanish n times in an interval 
in which the property W holds unless it is identically zero. 

For the vanishing of the coefficient 4; of Theorem II is precisely the 
condition that there exist a solution of (1) not identically zero passing 
through the » points (some of which may be coincident) involved in Aj. 


* A™-)(t) is not identically zero in any interval between 2, and 2y. 
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These » points were arbitrary in the interval so that the corollary is 
established. 

Now let us show that if A(¢) is a function of one sign in an interval 
a< x2 <b for every set of conditions (B) in this interval, then the property W 
holds in the interval a<2<b. We prove first that no solution of (1) 
can vanish m times in the interval a<a#<b. For suppose there were 
such a solution w. Let x, be the point of vanishing nearest b, and let «’ 
be a point between ay and 6. Now determine a solution « of the differ- 
ential system 

Lw = 1, 


w(x’) = wo = = (2’) = 
Form the function «(z) which is 


u(x) + Mw(z), 
u(x), 


where MM is a constant to be determined. We wish to show that M can be 
so determined that (a) vanishes n+ 1 times in <b. Now u(2’)+0 
since 2, <2’, and hence it follows that w(z’) has the sign of u(z’). 


Choose a point x” between x’ and b for which w(z”) is not zero. Such 
a point exists since w(x) = 0 in any interval. Now choose M so that 
u(x”) will have a sign opposite to that of w(a’); uw will then vanish 
between x’ and x”. wu has the required degree of continuity to apply the 
mean-value theorem, and A(¢) is a function of one sign. Hence Zu must 
change sign. But Zw is equal to zero in the interval a<x< 2’ and 
to M in the interval z’< x <b, and does not change sign. The contra- 
diction shows that w can not vanish » times in a<a2< b. 

But if no solution of (1) vanishes ” times in a < x<b, it is a simple 
matter to show that the property W holds in a<ca<b. For the principal 
solutions for the point a, 1%, us, --+, Un, are suitable functions. The Wrons- 
kian = W(m, Us, -++, ue) does not vanish in a<a<b. For suppose 
it vanished at a point ¢ of that interval. Then a function 


Cy Uy + Cota 


could be determined not identically zero and having k zeros in the point c. 
But this function would have » —k zeros in a, and a total of nm zeros in 
a =< «<b. This is impossible. We may now state the following 
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THEOREM III, A necessary and sufficient condition that the vanishing of 
a function u (with the required degree of continuity and not identically zero) 
at n+1 arbitrary points of an interval (a, b) should imply the change of 
sign of Lu at an intermediate point is that the property W hold in (a, b). 
A simple example will suffice to show that Theorem I is stronger than 
Theorem III. Take 
Lu = u"+u, 


= 0< 


sint 0 1 


A(t) = +sin(a,—?) sinx, cosx, 
| 


sina, COS | 
A(t) = sint sin(z,— 72), 


= sinz, sin(t—~ 7), 


Suppose now that and 
Then 
A(t)<0, a. 


A(t) is a function of one sign in the interval (0, x.) which may clearly be 
of length greater than a. (In fact it may be as near to 27 as we like.) 
Yet the property W can not hold in any interval of length greater than z, 
in as much as some solution of (1) will vanish twice in such an interval. 
This example suggests possible generalizations of Pélya’s results. 

It should be pointed out that Theorem I might easily be made to apply 
to the most general linear boundary conditions. 
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